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We classify distinct types of quantum number fractionalization occurring in two-dimensional topo- 
logically ordered phases, focusing in particular on phases with Z2 topological order, that is, on 
gapped Z2 spin liquids. We find that the fractionalization class of each anyon is an equivalence class 
of projective representations of the symmetry group, corresponding to elements of the cohomology 
group J/^(G, Z2). This result leads us to a symmetry classification of gapped Z2 spin liquids, such 
that two phases in different symmetry classes cannot be connected without breaking symmetry or 
crossing a phase transition. Symmetry classes are defined by specifying a fractionalization class for 
each type of anyon. The fusion rules of anyons play a crucial role in determining the symmetry 
classes. For translation and internal symmetry, braiding statistics plays no role, but can affect the 
classification when point group symmetries are present. For square lattice space group, time re- 
versal and S0(3) spin rotation symmetry, we find 2 098176 ~ 2^^ distinct symmetry classes. Our 
symmetry classification is not complete, as we exclude, by assumption, permutation of the different 
types of anyons by symmetry operations. We give an explicit construction of symmetry classes for 
square lattice space group symmetry in the toric code model. Via simple examples, we illustrate 
how information about fractionalization classes can in principle be obtained from the spectrum and 
quantum numbers of excited states. Moreover, the symmetry class can be partially determined 
from the quantum numbers of the four degenerate ground states on the torus. We also extend our 
results to arbitrary Abelian topological orders (limited, though, to translations and internal sym- 
metries), and compare our classification with the related projective symmetry group classification 
of parton mean-field theories. Our results provide a framework for understanding and probing the 
sharp distinctions among symmetric Z2 spin liquids, and are a first step toward a full classification 
of symmetric topologically ordered phases. 
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I. INTRODUCTION 

One of the characteristic features of topologically or- 
dered states of matter^^"^ in two dimensions is the pres- 
ence of anyons - quasiparticle excitations with non-trivial 
braiding statistics. Another important feature is quan- 
tum number fractionalization: if some degree of symme- 
try is present, the anyons can carry fractional quantum 
numbers. The charge e/3 quasiparticles of the v — l/i 
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Laughlin fractional quantum Hall state^ are a celebrated 
example of this phenomenon. The fractional charge of 
these excitations has been directly observed,^"'' while a 
direct, unambiguous measurement of their statistics re- 
mains elusive.^ As in this case, it is important to recog- 
nize that fractionalization may often be easier to detect 
than other characteristic features of topological order. 

Given the important role of fractionalization in topo- 
logically ordered states of matter, it is important to de- 
velop a better understanding of the interplay among sym- 
metry, fractionalization, and topological order. Many of 
the most basic questions along these lines are not well 
understood, for instance: Among states having the same 
topological order and the same symmetry, are there dis- 
tinct types of fractionalization that can be used to dis- 
tinguish phases? If so, how can we describe and classify 
distinct types of fractionalization? What types of frac- 
tionalization are consistent with a given type of topolog- 
ical order? In this paper, we answer these questions for 
one of the simplest types of topological order, namely 
the topological order of the deconfined phase of Z2 gauge 
theory,'' which we refer to as Z2 topological order. We 
will introduce the notion of fractionalization class of an 
anyon, which describes its characteristic type of fraction- 
alization. 

To be more specific, we shall confine our attention to 
two dimensions, to zero temperature, and to local bosonic 
models (i.e., spin models with finite-range or exponen- 
tially decaying interactions). For simplicity, we exclude 
the possibility of spontaneously broken symmetry. In this 
setting, states with Z2 topological order are referred to 
as gapped^" Z2 spin liquids. "'^^"^'^ Despite the restriction 
to two dimensions, it should be noted that nowhere will 
we assume a strict two-dimensional system. That is, our 
two-dimensional system may lie on the boundary of a 
gapped three-dimensional bulk. This point may have in- 
teresting implications for future work, and we return to 
it in Sec. IX. 

Our results can also be viewed through the lens of clas- 
sification of distinct phases of matter. More specifically, 
we may ask for a classification of all distinct phases of 
matter sharing a given fixed topological order and fixed 
symmetry group. In this situation. It is known that 
many distinct "symmetry enriched" topological phases 
exist. ^""^^ The distinctions among these phases disap- 
pear upon breaking of all symmetries, while the topolog- 
ical order is unaffected. 

We shall see that specifying fractionalization classes for 
each type of anyon defines a symmetry class, so named 
because it determines the action of symmetry on the 
topological degrees of freedom. Two states in different 
symmetry classes are distinct phases, and cannot be adi- 
abatically connected without closing a gap or breaking 
symmetry. This is only a partial classification of phases, 
though, because a given symmetry class may contain 
more than one distinct phase. Nonetheless, symmetry 
classification is a first step toward classification of all 
phases sharing a fixed topological order and symmetry 



group. We provide such a symmetry classification for 
gapped Z2 spin liquids, for an arbitrary symmetry group. 
It should be noted that the symmetry classification we 
give here is not complete; for simplicity, we do not con- 
sider cases where some symmetry operations permute the 
different types of anyons. This kind of action of the sym- 
metry group is "beyond fractionalization," and its study 
is left for future work. 

Our approach is not limited to Z2 topological order. 
Indeed, if only translation and internal symmetries are 
considered, we describe the straightforward extension of 
our results to arbitrary Abelian topological orders in 
Sec. V. Including full space group symmetry, with fur- 
ther work, we believe our approach could be extended to 
give a symmetry classification for an arbitrary Abelian 
topological order. We have not yet considered extensions 
to non- Abelian topological order or to topological order 
in three dimensions. 

Classifications such as ours are useful because they pro- 
vide a systematic basis for understanding sharp distinc- 
tions among phases. Along these lines, we hope that 
our results will be of use in finding new ways to identify 
and distinguish different spin liquids in numerical simu- 
lations and in experiments. Indeed, our focus on gapped 
Z2 spin liquids is partly motivated by the recent strik- 
ing evidence that such states are present in simple, fairly 
realistic S = 1/2 Heisenberg spin models on the J1-J2 
square lattice and the kagome lattice. '^""'^^ We do present 
some results here touching on determination of fraction- 
alization and symmetry classes in numerical simulations 
(Sees. IIIC and VII), but substantial further progress is 
likely possible, and we hope to stimulate further work in 
this direction. 

We begin with the familiar observation that quantum 
mechanics allows for symmetries to be realized projec- 
tively. One classic example is the fact that rotation by 
27r gives a phase —1 when acting on a wavefunction for 
a single half-odd integer spin. Another is the magnetic 
translation group of a single particle in a uniform mag- 
netic field, where two translation operators and Ty 
do not commute but instead satisfy T^Ty = e^'^TyT^- 
More generally, but somewhat loosely, we say symmetries 
are realized projectively when identities among group el- 
ements hold only up to a phase when acting on a quan- 
tum state. Group representations with this property are 
called projective representations. 

On the other hand, for any local bosonic model de- 
scribing a spin system built from an even number of elec- 
trons (or, for that matter, an even number of neutral 
atoms), symmetry operations act linearly — as opposed 
to projectively — on many-body wavefunctions. For in- 
stance, TxTy = TyTx- If one has a system with an odd 
number of electrons, we can always consider a larger sys- 
tem with an even number, so, with this constraint on 
lattice size in mind, we assert that symmetries act lin- 
early on the many-body wavefunctions of any physically 
reasonable local bosonic model. 

However, it is well known that, in general, symmetries 
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act projectively on anyons. For instance, S0(3) spin ro- 
tation symmetry acts projectively on the 5 = 1/2 spinon 
quasiparticles appearing in many gapped Z2 spin liquids. 
The crucial issue is how to describe and distinguish such 
projective actions to arrive at a set of fractionalization 
and symmetry classes. 

Before proceeding, we first have to briefly mention 
some facts about Z2 topological order. There are four 
particle types or classes of quasiparticle excitations, de- 
noted 1, e, TO and e. It may be helpful to think of these 
in terms of Z2 gauge theory coupled to bosonic matter 
fields; the deconfined phase of such a theory is a concrete 
realization of Z2 topological order. The e-particles are Z2 
electric charges, the m-particles are Z2 magnetic fluxes, 
and e-particles are e-m bound states. 1-particles, also 
referred to as "trivial" particles, are excitations that are 
not part of the topological structure. The non-trivial par- 
ticles {i.e., anyons) have non-trivial braiding statistics: 
any two distinct anyons {e.g. an e and an m) have 6 — tt 
mutual statistics. The e- and m-particles are bosons, but 
e-particles are fermions due to the mutual statistics of e 
and TO. 1-particles are bosonic, and have trivial mutual 
statistics with the other particle types. The fusion of two 
particles gives a unique third particle type, according to 
the fusion rules 

e X e = TO X TO = e X e = 1 (1) 
1x1 = 1, exl = e, TO xl = TO,exl = e, 
eXTO = e,exe = TO,TOXe = e. 

It is important to note that these properties are un- 
changed under the relabeling e <-> to. 

Now, to state our results, the action of the symmetry 
group on each type of topological quasiparticle is given 
by a projective representation, which is associated with 
a Z2 central extension of the symmetry group. (For 1- 
particles, this is always the trivial extension.) These cen- 
tral extensions can be grouped into equivalence classes, 
which we call fractionalization classes. Fractionalization 
classes are in one-to-one correspondence with elements 
of the cohomology group i?^(G', Z2). Symmetry classes 
are then defined by specifying the fractionalization class 
for each anyon. The symmetry class is a universal prop- 
erty of a Z2 spin liquid phase; that is, two states with 
different symmetry classes cannot be adiabatically con- 
nected without breaking symmetry. The fractionaliza- 
tion class for each anyon follows from the other two by 
fusion, so only two elements of iJ^(G', Z2) need be spec- 
ified. Equivalently, one can instead specify a single ele- 
ment of H'^{G,Z2 X Z2). 

Pairs of elements of H^{G, Z2) are not quite in one-to- 
one correspondence with distinct symmetry classes. This 
occurs because pairs of e and to fractionalization classes 
related by relabeling e -f-)- to are not distinct. 

If G consists only of translations and internal symme- 
tries, braiding statistics play no role in this classification. 
In this case, the fractionalization class of, say, e is given 
simply by the i/^(G, Z2) group product of the classes 
for e- and TO-particles. However, the statistics can enter 



when G contains more general space group operations, 
and in this case the 77^(G,Z2) product can be "twisted" 
by statistics. 

We hope that the reader is not discouraged at this 
stage by the appearance of perhaps unfamiliar mathe- 
matics. The necessary terminology and results are ex- 
plained in a self-contained fashion in Sec. IIIB. In our 
opinion, learning this material does not require any spe- 
cial mathematical sophistication. Group cohomology is 
certainly more sophisticated, but only the second coho- 
mology group appears, and that only as a convenient 
name for the group of equivalence classes of group exten- 
sions. 



A. Prior work 

The idea that symmetry acts projectively on topolog- 
ical quasiparticles also lies at the heart of X.-G. Wen's 
projective symmetry group (PSG) classification of mean- 
field spin liquid states, which is a key inspiration for 
our work and can be viewed as an attempt to answer 
some of the same questions. However, PSG classifica- 
tion, while a very useful tool, does not give a symmetry 
classification. PSG classification begins with a parton 
construction, where, for instance, the spin operator is 
written as a bilinear of bosonic or fermionic spinon oper- 
ators. One then constructs a mean-field theory in terms 
of the partons, and such mean-field theories are classi- 
fied by PSG. Fluctuations about mean-field theory can 
be incorporated by coupling the partons to a dynamical 
gauge field, giving a true low-energy effective theory. 

The PSG classification is inherently tied to parton ef- 
fective theories. A symmetry classification should be 
built only on the essential defining properties of Z2 topo- 
logical order — namely, the types of anyons and their fu- 
sion and braiding properties. Parton constructions pro- 
vide a concrete means to realize these properties, but 
there is no reason to believe they do more than this. Put 
another way, Z2 topological order does not seem to be 
essentially linked to parton theory, so in our view par- 
ton theory and PSG do not provide the right language to 
construct a symmetry classification. We provide a more 
detailed discussion contrasting PSG classification with 
our symmetry classification in Sec. VIII. 

We also note that some ideas having significant over- 
lap with ours were outlined previously by A. Kitaev, in 
Appendix F of Ref. 36. In particular, taking the liberty 
of translating results presented there into the language 
of this paper, it was asserted that for a general topolog- 
ical order the symmetry classes are given by elements of 
-ff^(G, r2), where G is the symmetry group and r2 is a 
finite Abelian group determined by the type of topologi- 
cal order. For Abelian topological orders, r2 is the group 
of fusion rules, "^^ so r2 = Z2 x Z2 for Z2 topological or- 
der, agreeing with our results. In fact, for an arbitrary 
Abelian topological order, we show in Sec. V that our 
approach also reproduces Kitaev's result if G consists 
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only of translations and internal symmetries. This may 
also hold for more general space group symmetries, but 
there are subtleties having to do with the role of braiding 
statistics that we have only addressed for Z2 topological 
order. 

A number of other prior works have also investigated 
related questions. ^^"^^'^^''^'^ In particular, in Ref. 21, the 
idea of using a pair of PSGs, one for Z2 charges and 
one for Z2 fluxes, was introduced. This idea enters our 
symmetry classification via the need to specify two frac- 
tionalization classes (for instance, the e and m fraction- 
alization classes). Ref. 21 also showed that distinct pairs 
of charge and flux PSGs can be realized in the toric code 
model, in close connection to our analysis of the same 
model in Sec. VI. 



B. Outline 

We begin in Sec. II A with a review of Z2 topological 
order in two dimensions, introducing many of the basic 
ideas important for our symmetry classification, as well 
as much of the notation used in the rest of the paper. 
Of particular importance is the concept of superselection 
sectors. Next, in Sec. II B we briefly review the toric code 
model, the simplest concrete realization of Z2 topolog- 
ical order. 

Sections III and IV present the central results of the 
paper. In Sec. Ill A, the notion of fractionalization class 
of an anyon is introduced, focusing on the case of trans- 
lation and internal symmetry. The notions of symme- 
try localization and one-particle symmetry operators are 
also introduced, and play a crucial role. (One techni- 
cal detail is relegated to Appendix A.) We show that 
for translation symmetry alone there are two fraction- 
alization classes, and similarly for U(l) or S0(3) sym- 
metry. Next, in Sec. IIIB we introduce the mathemati- 
cal language needed to describe fractionalization classes, 
followed by a general discussion of the structure of frac- 
tionalization classes in Sec. Ill C. In Sec. Ill C, via simple 
examples, we also explain how fractionalization class in- 
formation can manifest itself physically in the excitation 
spectrum and quantum numbers of excited states. As 
part of this discussion, we introduce a "coarsened" clas- 
sification by Ut(1) fractionalization classes, which reflect 
information that is in a sense physically simpler than that 
contained in the full classification. Because point group 
operations can move anyons large distances, full space 
group symmetry requires the further considerations of 
Sec. HID. Finally, in Sec. HIE we work out the fraction- 
alization classes for the example of square lattice space 
group, time reversal and spin rotation symmetry, show- 
ing there are 2^^ such classes. 

In Sec. IV A, we describe our symmetry classification 
of Z2 spin liquids, which amounts to specifying fraction- 
alization classes for the e and m anyons. The crucial issue 
is to determine how the e fractionalization class follows 
from the e and m classes. Following a discussion of the 



counting of distinct symmetry classes, we move on to the 
case of translation and internal symmetry in Sec. IV B, 
where we show that the e class is given simply by the 
i?^(G, Z2) group product of the e and m classes. We 
also describe the symmetry classes for the case of trans- 
lation symmetry alone, and for SO (3) spin rotation alone. 
Section IV C discusses symmetry classes for space group 
symmetry, where the mutual statistics of e and m parti- 
cles leads, in general, to a twisting of the group product 
determining the e fractionalization class in terms of the 
e and m classes. We explicitly work out this twisting for 
the square lattice space group. 

In Sec. V, we extend our results to general Abelian 
topological orders for the case of translation and inter- 
nal symmetry. In Sec. VI, we explicitly construct one- 
particle symmetry operators for the generators of the 
square lattice space group in the toric code model, and 
show that three symmetry classes can be realized there, 
by tuning the signs of the two terms in the Hamiltonian. 
Section VII shows that some of the symmetry class in- 
formation can be extracted directly from the quantum 
numbers of degenerate ground states, as illustrated for 
the case of translation symmetry alone. We conclude 
with a comparison beween our classification and PSG 
classification in Sec. VIII, and a discussion of open issues 
and future directions in Sec. IX. 



II. REVIEW: Z2 TOPOLOGICAL ORDER 
A. General discussion 

Here, we review Z2 topological order. We employ the 
language of topological quasiparticle types and the asso- 
ciated superselection sectors. '^^ The notion of topological 
superselection sectors is particularly important for our 
classification. We begin by introducing these notions ab- 
stractly, and then, in Sec. II B, illustrate them using the 
exactly solvable Kitaev toric code model. The discus- 
sion here focuses on Z2 topological order, but can be 
generalized to arbitrary Abelian topological orders. 

We are concerned with local bosonic lattice models 
with an energy gap in two dimensions. One way to char- 
acterize Z2 topological order is by properties of excita- 
tions above the ground state. As discussed briefly in 
Sec. I, there are four topological particle types, denoted 
1, e, m and e. The e, m and e particles are anyons, while 
the 1-particles ("trivial" particles) are not. Under ex- 
change, all the particle types are bosons except e, which 
is a fermion. Excluding 1-particles, any pair of distinct 
particles have 9 — -k mutual statistics. 1-particles have 
trivial mutual statistics with the other particle types. 

The particles also obey the fusion rules given in Eq. (1) 
For example, the e x e = 1 fusion rule expresses the fact 
that two nearby e-particles can be viewed as a single 
1-particle. Because the fusion and braiding rules are in- 
variant under e o to, we are always free to relabel e ^ m 
if we wish. Other Abelian topological orders can be de- 





FIG. 1. (a) Two e-string operators (solid lines) with a single 
crossing point commute, (b) e-string (solid line) and m-string 
(dashed line) operators with a single crossing point anticom- 
mute. 



FIG. 2. (a) An e-string (solid line) and m-string (dashed 
line) can be viewed in combination as an e-string. The arrow 
points from the m-string toward the e-string. (b) Twisted e- 
string where the e-string passes underneath the m-string. (c) 
Twisted e-string where the e-string passes over the m-string. 
Note that the configurations in (b) and (c) are related by a 
minus sign. 



scribed in the same way; that is, one specifies a set of 
particle types, fusion rules, and both exchange and mu- 
tual statistics. 

At this point, it is useful to introduce some terminol- 
ogy. We consider a region R defined as some subset of 
all lattice sites. Without worrying too much about preci- 
sion, we also assume that R is has no small holes or rough 
edges. That is, we want to be able to define R by draw- 
ing one or more sufficiently smooth boundary curves, and 
selecting the lattice sites in the interior. We will almost 
always assume i? is a union of disjoint simply connected 
regions. We denote the complement of R by R. The full 
Hilbert space is the tensor product ■H(i?) (8) 'H{R), where 
'H{R) is the Hilbert space of region R. We say an opera- 
tor O is supported on R if it can be written O = Or^^r- 
That is, if O is supported on R it may act non-trivially 
on R, but acts as the identity operator on R. 

It is a crucially important defining property that no lo- 
cal operator can create a single isolated anyon. However, 
for instance a pair of e-particles can be created locally 
due to the fusion rule e x e = 1, since isolated 1-particles 
can be created locally. Two e-particles created in this 
way can then be separated to obtain isolated e-particles. 
This separation can be accomplished by acting with a 
string operator, which is supported on a linear region 
connecting the initial and final positions of a single e- 
particle. String operators only modify locally observable 
properties of a state on which they act near the ends 
of the string. That is, there is no way to discern that a 
string operator has been applied to some state by making 
local measurements along the length of the string, away 
from the ends. String operators need not have ends, and 
can form closed loops. Noncontractible closed strings are 
related to the topological ground state degeneracy, which 
we discuss below. 

There are three distinct types of string operators, asso- 
ciated with the three types of anyons. For instance, the 
string operators associated with e-particles are referred 
to as e-strings. Two string operators of the same type 
commute (Fig. la). This holds as long as their ends are 
well separated; if that is not the case, the commutation 
relations will depend on details of the ends. On the other 
hand, two string operators of different types anticommute 



if they cross an odd number of times and commute if they 
cross an even number of times (Fig. lb). This property, 
which follows from the mutual statistics, can be expressed 
by saying that we get a minus sign whenever we move a 
string of one type through a string of another type, at a 
single crossing point. In Fig. 1 and throughout the pa- 
per, we adopt the graphical convention that strings are 
drawn on top when the corresponding operator lies to the 
left in a product of operators. That is, the strings on the 
bottom act first on a state, followed by strings on top. 

An e-string can be viewed as a pair of nearby, parallel 
e and m strings, as depicted in Fig. 2a. Since the e-string 
has to lie on one side or the other, e-strings thus carry 
an orientation (Fig. 2a). The orientation changes at a 
point where one of the e or m constituent strings passes 
under the other; at such a point we say the e-string is 
twisted. There are two kinds of twists, since the e-string 
can pass over or under the m-string (Figs. 2b, 2c); these 
two twists are related by a minus sign. 

We will primarily be interested in considering Z2 topo- 
logically ordered states on a torus {i.e. with periodic 
boundary conditions). In this case, we can form non- 
contractible loops in both x- and y-directions. Let C% be 
a closed e-string operator winding once around the sys- 
tem in the cc-direction, with corresponding definitions for 
Cy, and £™. We can think of the product C%C™ as a 
e-string running in the a;-direction, so it is not necessary 
to introduce more operators to represent non-contractible 
e-strings. These string operators satisfy the following 
commutation and anticommutation relations: 

{C^, Cy } {-Cj; ^ 0, 

\C^, Cy] [C^ , Cy ] [C^, C^ ] — [Cy, Cy ] 0. (2) 

We also assume 

iClf = {Clf = {C7f = {C^r = 1- (3) 

This can be justified by noting that, for instance, mov- 
ing a pair of nearby e-particles (equivalent to a 1-particle) 
around a closed loop should be equivalent to doing noth- 
ing, except perhaps accumulating a phase that can be 
removed by a trivial redefinition of the string operators. 



We refer to Eqs. (2) and (3) as the loop algebra. This 
algebra has a single, four-dimensional irreducible repre- 
sentation, and this implies that the ground state on a 
torus must be fourfold degenerate. We note that there 
are situations where the loop algebra and the topologi- 
cal ground state degeneracy are both modified by choice 
of boundary conditions,'*" but we will not consider such 
cases. 

Associated to each of the four particle types is a topo- 
logical superselection sector. To understand what this 
means, it is first helpful to think about a single isolated 
and localized e-particle in an infinite plane, where no 
other excitations are present. Starting from such a state, 
we define an arbitrarily large but finite connected region 
R*^ containing the e-particle. It is important that the 
system be "locally in the ground state" near the bound- 
ary of i?*^ (and outside of i?*^), in the sense that local 
measurements in these areas should give the same result 
as in the ground state. Then we can obtain all states 
in the e-sector by acting with (almost)'** arbitrary oper- 
ators supported on R^. The resulting states correspond 
to moving the e-particle to different positions, "dressing" 
it in various ways, modifying any internal quantum num- 
bers it may carry, and so on. The e-sector is thus closed 
under action of operators supported on i?*^. Moreover, 
no operator supported on can act on an e-sector state 
and turn it into a state belonging to a different superse- 
lection sector. 

We will also have occasion to consider regions that are 
not connected. To handle this situation, we make the 
following definition: a s-operator on R is an operator 
that, restricted to R, consists entirely of string operators 
(of any type) connecting disconnected components of R. 
If an s-operator on R has only, say, e-strings in R, we 
call such an operator an e-operator on R. Again, if the 
region R'^ contains an isolated e-particle and no other ex- 
citations, we can obtain all e-sector states (for the region 
R*^) by acting with s-operators on i?*^, and the e-sector 
is closed under the action of such operators. 

This discussion needs to be modified on a finite torus, 
where non-trivial particles must occur in pairs. Indeed, 
in this situation all physical states belong to the 1-sector, 
because any state can be obtained from a ground state 
by acting with operators supported on the whole system. 
To see that the superselection sectors still have mean- 
ing here, consider a state with two localized and well- 
separated e-particles, with no other excitations present 
(Fig. 3a). Any operator acting on a ground state to cre- 
ate such a two-particle state will involve an e-string con- 
necting the positions of the two particles; therefore, it is 
useful to think of the particles as being connected by an 
e-string. 

Now, as illustrated in Fig. 3a, we define two regions R\ 
and i?| by drawing a box around each e-particle. The 
boundaries of these regions, and the space outside the 
regions, should be locally in the ground state. From such 
a reference state, acting with arbitrary operators sup- 
ported on, say i?f , one generates all e-sector states in the 
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FIG. 3. (a) State with two isolated e-particles, with e-sector 
regions R1 and R^. The e-string connecting the particles is 
shown as a solid line. The two regions can be combined to- 
gether to give the 1-sector region . (b) State as in (a), but 
where the e-particle in _Ri has split into isolated m and e par- 
ticles as allowed by the e = m x e fusion rule. R\ can thus 
be subdivided into _R™ a^nd R\ as shown. Strings connecting 
the particles are not shown. 



FIG. 4. Thick bonds depict the four edges meeting the vertex 
s and bounding the plaquette p. 



region i?^. 

More generally, we can decompose a state into regions 
R], Rf, RY^, and Rj. Such a decomposition is not unique 
and can be modified according to the fusion rules. For 
instance, going back to the example of two e-particles, 
if we draw a larger box containing both R^ and R2, the 
resulting new region i?* contains a state in the 1-sector 
due to the fusion rule e x e = 1. On the other hand, due 
to the e — m X e fusion rule, there can be a state in R^ 
consisting of a localized m-particle well-separated from 
a localized e-particle, with no other excitations. In this 
case we can subdivide Rf into and i?™, as shown in 
Fig. 3b. 
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FIG. 5. Depiction of contours (thick solid line) and C™ 
(thick dashed line) used to define the loop operators C% and 
£™ , respectively. 



B. Toric code model 

The toric code is a spin model that makes manifest all 
the essential features of the Z2 theory with a minimum 
of extra structure. It is therefore a very useful testbed 
for our discussion, and in section VI we will present ex- 
plicit constructions that give substance to the general 
considerations of Sections III and IV. The model can be 
defined on any lattice in two dimensions, but we restrict 
ourselves to the square lattice. 

The model consists of spin- 1/2 degrees of freedom on 
the edges of the square lattice. We label lattice sites by r, 
and write Pauli matrices acting on the nearest-neighbor 
edge (r, r') as cr^r'' ^^'^ Four edges meet at each 

site to form a vertex, denoted s, and four edges bound 
each plaquette, denoted p (Fig. 4). The Hamiltonian is 
built from the following products, associated respectively 
with vertices and plaquettes: 

n <r', B,= II (4) 

(r,r')gs (r,r')ep 

These operators can be viewed as measuring Z2 charge 
and flux, respectively. The Hamiltonian is 



the X, y directions, where the C™^ contours pass through 
perpendicular edges. It is straightforward to check these 
operators satisfy the loop algebra. More generally, e (m) 
strings are products of af (af ) along an appropriate con- 
tour. 

We can be even more explicit by constructing the state 



1^0) = [n 



1 



(1 + ^) |{af = 1}). 



(8) 



Htc = ~Ke '^As- K,n Ap. 



(5) 



This is easily seen to be a ground state, and £^^y\tpo) = 
\ipo}- However, this state is not an eigenstate of -C^y, 
and the other three ground states are £™|V'o)j 'C™|V'o) 

and r^r^iVo). 

Along the same lines, we can also construct con- 
tractible e- and m-strings. Any such e-string can be writ- 
ten as a product of Bp operators; a single Bp operator 
is an elementary contractible e-string encircling a single 
plaquette. Similarly, contractible m-string operators are 
products of As operators. Therefore, the ground states 
are eigenstates of all contractible string operators, with 
eigenvalue unity. 

There is a gap to excited states, which have vertices 
where Ag — —1 and/or plaquettes where Bp — — 1. 
Stars with Ag = —1 are e-particles, and plaquettes with 
Bp = —1 are m-particles. (Remember, that which we 
call e and which m is arbitrary.) To create, for instance, 
a state with two isolated e-particles, one can act on a 
ground state with a product of erf Pauli matrices (an 
e-string) on a contour connecting the desired particle po- 
sitions. Based on the above discussion, it is clear that we 
can "slide around" the string connecting the two parti- 
cles with no effect whatsoever on the state. This state- 
ment needs to be weakened slightly when we consider 
states with both e and m-particles; in such cases moving 
strings around can change the state by a minus sign, ei- 
ther when we bring an e-string through an m-string, or 
when a string of one type "slides over" a particle of the 
other type. Still, the string positions are clearly unob- 
servable. 



III. FRACTIONALIZATION CLASSES 



The exact eigenstates of Htc are easily constructed, be- 
cause [Htc, As] = [Htc, Bp] = [As, Bp] = 0. 

Here we assume Ke,K„i > 0, although in Sec. VI we 
consider more general situations. Any ground state sat- 
isfies As = Bp = 1. On a finite torus there are four such 
ground states, which can be seen by explicitly construct- 
ing the loop operators as products of Pauli matrices. 



n 



x,y J_ J_ r,r' ' 

(r,r')J-C- 



(6) 
(7) 



as illustrated in Fig. 5. Here the contours C| y consist 
of a path of lattice edges winding around the system in 



In this section we will discuss and classify the action 
of symmetries on a single type of anyon in the Z2 the- 
ory. Most of the discussion will focus on e-particles, with 
occasional comments on e-particles, due to the different 
nature of e-strings. Everything we say also clearly holds 
for TO-particles, since the topological properties do not 
change under relabeling e O m. 

We assume that symmetry operations do not change 
one type of anyon into another. The symmetry group 
may consist of internal symmetries (including anti- 
unitary time reversal), translation symmetry, and general 
space group operations. We begin by introducing the no- 
tion of fractionalization classes for translation and inter- 
nal symmetry (Sec. HI A). Next, introducing the mathe- 
matics of group extensions and their equivalence classes 
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FIG. 6. Depiction of the state \ipa) with two e-particles. The 
operator Oa creating this state is supported on the two shaded 
circular regions and along the solid line connecting them. Oa 
is an e-string operator along the solid line, which is referred 
to as LI outside of the regions _Rf . Oa is thus an e-operator 
on the union of the R1. 



(Sec. IIIB), we discuss the general structure of fraction- 
alization classes (Sec. Ill C). We then show that the same 
general structure continues to hold for space group sym- 
metry in Sec. HID. 



A. Translation and internal symmetry 

To introduce the notion of fractionalization classes, we 
begin with translation symmetry, then argue that the 
same structure holds for internal symmetry (including 
time reversal). We discuss e-particles for concreteness, 
but all statements apply just as well to e-particles, except 
where explicitly noted. 

We consider translation symmetry generated by and 
Ty, satisfying the relation T^TyT'^T'^ — 1, which holds 
acting on all physical states (see Sec. I), and thus on 
the 1-sector. We wish to understand the action of trans- 
lations on states with two localized, well separated e- 
particles. More formally, we consider a family of states 
{IV'a)} (labeled by a), that can be decomposed into two 
fixed, connected e-sector regions [i = 1,2), as de- 
scribed in Sec. II A, and are otherwise locally in the 
ground state (Fig. 6). We also assume that |V'q) — 
CalV'o), where j^/'o) is a ground state. The operator Oa is 
an e-operator on the union of R\ and i?2 , with an e-string 
running between the two components. For simplicity, we 
assume that |?/'o) satisfies T!j:|-0o) = Ty\iljQ) = iV'o)- This 
assumption is not necessary and we describe how it can 
be relaxed in Appendix A. Any desired combination of 
e-sector states in the two regions Rl can be produced by 
the above construction. 

We proceed by making the crucial assumption that 
symmetry operations can be localized to regions sur- 
rounding the excitations in the states We refer to 
this property as symmetry localization. The basic idea is 
that a symmetry operation (such as translation) changes 
the e-particle state in one of the regions to another e- 
particle state in the same region, and that it should be 
possible to accomplish such a change locally. Formally, 

TM=n{^)Tm\^o.). (9) 



where T^{i) is supported on R^. The operators T|(i) are 
independent of a. The corresponding statements are also 
assumed for Ty. The operator T!|(l), for instance, can 
be interpreted as a "one-particle" translation operator, 
that translates the e-particle in region i?f against the 
translation-invariant medium of the ground state |V'o)- 
It is straightforward to generalize this discussion to a 
state with multiple e-sector regions. We further examine 
and justify the assumption of symmetry localization at 
the end of this section. 

We note that symmetry localization fails when a sym- 
metry operation changes one type of anyon into another 
type, because this cannot be accomplished by any local 
operator. However, we exclude this situation by assump- 
tion. 

Now, we have 

= n (T^.^«r;(z)K(z)]-Mr;(*)]-^)|^a). 

1=1,2 

For this to hold for all a, wc must have 

T^A^K{i)[T^Ai)r'[T^M-' - e'^^ (11) 

If this were not true, there would be a state Ii/jq) on 
which the identity operator acts nontrivially. We have 
gi0ig«02 = \^ which is a consequence of the fusion rule 
e X e = 1. Since there is no difference between the 
two regions, we also expect e*"^^ = e*'*^. To show this, 
consider instead a state with four e-sector regions i?f 
[i = 1, ... ,4), with e"^' defined as above. Then for any 
pair i ^ we can fuse Rl and i?| to obtain a 1-sector 
region, acting on which we must have TxTyT~^T~^ = 1, 
implying e^'^'e*'^^ = 1. This implies the e''^' are all equal 
and e*"^' — ±1. Therefore, dropping the label distinguish- 
ing the two regions, we write 

T^^T^yinr^T^y)-^ ^ cjl^.y = ±1. (12) 

The parameter a^^^y defines the fractionalization class 
of the e-sector. Evidently, there are two fractionalization 
classes in the case of translation symmetry alone. It is 
important to emphasize, as follows from the discussion 
above, that a^'i^^y is constant on the e-sector. Putting es- 
sentially the same argument in different terms, suppose 
there is one type of e-particle with cr^^^^y = 1 and another 
with al^iy — —1. We could then fuse these to obtain 
a 1-particle acting on which T^TyT^^Ty^ = —1, a con- 
tradiction. Since crfj-tj, is discrete and is constant on the 
e-sector, it cannot change within a Z2 spin liquid phase, 
so long as translation symmetry is preserved. Therefore, 
crfjjjy is a universal property of Z2 spin liquids with trans- 
lation symmetry. 

There is some arbitrariness in the definition of T!j.. 
Looking at Eq. (9), clearly we can redefine T|(i) — >■ 
—T^{i) with no physical effect. More generally, the re- 
definition T!^{i) e"t'T^{i) sends e^-^T^, acting 
on a state with n e-particles. This transformation should 
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leave unchanged (apart from possible overall multipli- 
cation by a phase, independent of n), which only occurs 
when (j) — 0,77, in which case T^. — ^ Tj. since n is even. 
Therefore we are allowed to redefine T^{i) — > —T^{i), and 
similarly Ty{i) — >■ —Ty{i). Note that such redefinitions 
do not affect (Jtxty 

Now wc generalize the above discussion to the case of 
unitary internal symmetry, perhaps also combined with 
translation symmetry. By internal symmetry, we roughly 
mean any symmetry operation that does not move the 
lattice. This includes on-site symmetries such as spin 
rotation, but we need not limit ourselves to strictly on- 
site symmetry. If 5 is a symmetry operation, we again 
assume symmetry localization, that is 

^|Va)=^'=(l)^'^(2)|V'a), (13) 

where S'^{i) is supported on R^. The logic is identical to 
the case of translation: it should be possible to accom- 
plish the operation S by making local modifications near 
the two quasiparticle excitations. Again, we are free to 
redefine S"^ — ?► —S"^ with no effect on the physics. Sup- 
pose we have a relation among symmetry operations of 
the form S1S2 ■ ■ ■ Sk = 1- Then, following the arguments 
above, 

S1S^,---S^,^±1. (14) 

Specifying such Z2-valued parameters for all group rela- 
tions among symmetry operations specifies the fraction- 
alization class of the e-sector. At the present stage of 
the discussion, it may not be clear how to make this 
last statement precise. This can be accomplished in a 
straightforward fashion after the discussion of the follow- 
ing section, where the mathematics of group extensions 
and their equivalence classes is introduced. 

It is worth explicitly discussing the particularly simple 
and familiar cases of U(l) and S0(3) internal symme- 
try. In the case G = U(l), 1-sector states {i.e., physi- 
cal states) carry integer U(l) charge, or more generally 
are superpositions of states with different integer U(l) 
charges. Alternatively, denoting with R{0) a U(l) ro- 
tation by angle 6, we have R{2-k) = 1. On the e-sector, 
however, we may have i?*^(27r) = ±1, corresponding to in- 
teger (-1-1) and half-odd integer (—1) U(l) charge. These 
are the only two fractionalization classes. For instance, 
e-particles cannot have other charges {e.g., 1/3 charge), 
since combining two e-particles must always give an in- 
teger charge due to the e x e = 1 fusion rule. Moreover, 
e-particles with charge 1/2 and charge 3/2 are not dis- 
tinct classes; starting with a charge-1/2 e-particle, one 
can fuse it with a charge-1 1-particle to obtain a charge- 
3/2 e-particle. Therefore, charge-1/2 and charge-3/2 e- 
particles always appear together in the spectrum. This 
example points out that fractionalization classes are not 
simply distinct irreducible representations of the symme- 
try group, and are instead a coarser type of classification. 
The situation for G — SO (3) spin rotation is similar. De- 
noting with Rs{9n) spin rotation by 6 about the n-axis. 



on the e-sector we have Rl{2iTn) = ±1, corresponding to 
the two fractionalization classes of integer spin (-1-1) and 
half-odd-integer spin (—1). 

Finally, we consider the case of anti-unitary time re- 
versal, which can be written T = UtK, where K is com- 
plex conjugation and Ut is a unitary operator. Complex 
conjugation is a global operation on a state and cannot 
sensibly be localized to a region, so in this case symmetry 
localization takes the form 

T\iJa) = U^{l)U^{2)K\ij^), (15) 

where C/t(*) is supported on Rf. The relation = 1 
implies 

U^{U^y EE {T^' f = ±1. (16) 

Here, in the interest of concise notation, we have made a 
formal definition of {T^)^. 

We can also consider relations involving time reversal 
and unitary symmetry operations. For instance, suppose 
TST~^S~^ = 1 for some symmetry operation S. This 
implies 

u^{sy{u})-\sT^ = t^s^t'THst^ = ±i, 

(17) 

where again the expression involving T*^ is a formal def- 
inition. 

We now discuss the assumption of symmetry localiza- 
tion in more detail. Let R be the union of the i?| regions, 
and R the complement of R. In addition to acting on the 
e-particles in R with some symmetry operation S, we also 
have to act on the e-string operator in R. Now, moving 
or otherwise modifying the string can only result in an 
overall phase, but the question is then whether this phase 
can always be factorized into a product of two phases as- 
sociated with each region. In general, we expect 

S\i;^) = S^{l)S^{2)L'{l,2)\i;^), (18) 

where L'^{1, 2) is a closed e-string loop that accomplishes 
the necessary transformation of the e-string in R (see 
Fig. 7a). Restricted to R, L^{1,2) = L)L1, where if is 
identical to Oa restricted to R; that is, LI is the e-string 
of the "initial" state \ipa)- is the desired "final" state 
e-string. Since L''{\, 2) is a contractible loop operator, it 
can be written as a product of smaller contractible loops 
of e-string as shown in Fig. 7b. We assume that e-string 
operators square to unity, so the elementary contractible 
loop operators have eigenvalues ±1. Assuming some de- 
gree of spatial homogeneity, by combining and splitting 
the elementary loops as needed, it should be possible to 
choose all the elementary loops to have eigenvalue unity, 
except possibly those at the ends. Therefore, L^(l,2) 
can be broken in the middle and deformed as shown in 
Fig. 7c, without accumulating any phase factors. The 
remaining loops on the two ends can then be absorbed 
into S"^(l) and >5"^(2), and symmetry localization holds. 

The discussion of transforming the string also goes over 
to the case of e-particles, where the string is oriented. It 
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(b) 












(c) 




FIG. 7. Depiction of the action of a symmetry operation S on 
a two e-particle state, where S is either translation or an in- 
ternal symmetry operation. The initial state on which 5* acts 
is shown in Fig. 6. S acts nontrivially in the shaded circular 
regions, and also as the closed e-string loop (solid line), as 
shown in (a) . Outside the regions Rl , Lf is the e-string of the 
initial state \tpa}, and LJ is the e-string of the transformed 
state S\'ijja}- In (b), the e-string loop is divided into a prod- 
uct of smaller loops, chosen so that each gives unity acting on 
the ground state, except possibly the two loops at the ends. 
Therefore, the loops in the center can be eliminated (c), leav- 
ing only two loops contained entirely within the regions R^, 
which can be absorbed into the definition of 5"^(i). 



is still possible to construct the necessary closed e-string 
loop operator L'^(l,2) as above for e-particles. The only 
difference from the case of e-particles is that, depending 
on the relative orientations of the and strings, it 
may be necessary for the i^(l, 2) to be twisted near the 
ends, inside the two e-regions. Otherwise, the discussion 
proceeds exactly in the e-particle case. 



B. Group extensions and their equivalence classes 

Here, we give an account of group extensions and their 
equivalence classes. For the most part, we find it clearer 
to separate the mathematics from the physics, so the 
discussion here focuses on the mathematics. In learning 
this mathematics, we found it useful to consult Refs. 42- 
47. Application of the mathematics to the physics of 
fractionalization follows in Sec. IIIC. 

Consider a group G with elements g ^ G. We consider 



a projective representation, where the group element g is 
represented by a unitary matrix T{g). (We include anti- 
unitary group operations below.) When multiplying two 
F's we have 



r(.gi)r(.g2) = w(5i,ff2)r(gig2), 



(19) 



where ^(51,52) G U(l) is a phase. The presence of 
these U(l) phases is what it means for the representation 
to be projective. An "ordinary" representation where 
'^{91 1 92) = 1 for all (71,92 is referred to as a linear repre- 
sentation. We allow for the possibility that the projective 
representation F may be a linear representation; that is, 
any linear representation is a projective representation, 
but not vice-versa. To connect with the discussion of 
Sec. Ill A, F arises physically as the action of the sym- 
metry group on one of the superselection sectors. 

We restrict uj{gi,g2) G A, where ^ is a subgroup of 
U(l). In the physical applications of this paper, we will 
be interested in the case A = 1^2- The function a; (171, 52) 
satisfies an associativity constraint, because 



F(gi)F(52)r(53) 



^(.9i,52)'^(5i32,ff3)r(5ig253) (20) 
'^(Si, 5233)^(32, ff3)r(5ig253), 



where the two results are obtained by the two different 
ways of using associativity to evaluate the product of 
three F's. The associativity constraint is then 



'^(31,52)^(3152,. 93) = t^(5i, 5253)^(52, 53)- 



(21) 



Any function 0^(91,52) € A satisfying the associativity 
constraint is called a factor set, or sometimes an A-factor 
set. 

If G includes anti-unitary operations, and if A ^ Z2, 
the above discussion needs to be modified. If g is anti- 
unitary, then so is Y{g), which acts nontrivially on ele- 
ments of A by complex conjugation. For example, 

r(5i) w(5r2,53)r(5253) = ^"^(52,53)r(5i)r(52.g3) 

(22) 

for gi anti-unitary. This modifies the associativity con- 
straint. We define s{g) = 1 for g unitary and s{g) = — 1 
for g anti-unitary, and then 

w(5i,52)w(5i52,53) ^ ^{gi.g2g'i)i^''^^^\g2,gz)- (23) 

To remind ourselves of the nontrivial action of anti- 
unitary operations, we will call such factor sets Ar-factor 
sets. Since we are mostly interested in yl = Z2, this com- 
plication is not relevant to much of our discussion. How- 
ever, we will also have occasion to consider UT(l)-factor 
sets, so we will include the possibility of nontrivial action 
of anti-unitary operations in the discussion below. 

If 1^0(51,52) and 0Jb{gi,g2) are factor sets, then 
Wah(5i,32) = LOa{gi,g2)u}bigi,g2) is also a factor set. It 
is simple to check that this product defines an Abelian 
group structure on the collection of all factor sets. The 
product of factor sets is associated with tensor prod- 
ucts of representations: if Wa,Wf, are the factor sets of 
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the representations FajTb, respectively, then the prod- 
uct ujab = WaWb is the factor set of the tensor product 
representation (81 Ff,. 

Suppose that we allow a redefinition of the F's by 

r'ig) = A(.9)F(g), (24) 

where X{g) £ A. This induces the following transforma- 
tion of the factor set: 

^'(51,52) = A(.gi)A^(»^)(.92)A(.9i.g2)-ic.(5i,52). (25) 

oj' is also a factor set {i.e., satisfies the associativity con- 
straint). Two factor sets to and uj' are said to be equiv- 
alent if they are related by Eq. (25) for some X{g), and 
in this case we write w ~ w'. This notion of equivalence 
is reflexive (w ~ w), symmetric (uj' ~ w if a; ~ w') and 
transitive (if lo ^ u' and uj' ^ uj" , then lj ^ w"), there- 
fore ~ defines an equivalence relation that partitions the 
set of factor sets into equivalence classes. We denote the 
equivalence class of u by c{uj). Note that c(aj) = c(a;') 
if and only if cj ~ w'. Given a class c{uj), we say w is a 
representative of the class. 

The equivalence classes themselves form an Abelian 
group, with product defined by 

C{u}i)c{uj2) = c{uJiUJ2)- (26) 

This product is well-defined, in the sense that it does not 
depend on the representatives we choose for each class. 
The Abelian group of factor set equivalence classes is 
isomorphic to the cohomology group H^{G, At). If anti- 
unitary operations are not present, or if they act trivially 
on A as when A — Z2, we leave off the T subscript and 
write H^{G,A). We shall not bother to give a definition 
of H^{G, At) in terms of group cohomology, because, for 
our purposes, it is sufficient to view the group of factor 
set equivalence classes as the definition of H'^{G, At) (see 
footnote 17 of Ref. 48). We shall often refer to factor set 
equivalence classes as cohomology classes [i.e., elements 
of H^iG,AT)]. 

With this discussion behind us, we note that the pro- 
jective representation F (with anti-unitary operations) 
is associated with an Ay-extension of the group G. 
Roughly, such an extension is a new group E in which it 
makes sense to multiply elements of A and elements of 
G. The advantage of a defining an extension is that it is 
characterized entirely by G, A and uj, so we can equiva- 
lently speak of classifying factor sets or classifying group 
extensions. Formally, an Ar-extension is a group E such 
that: (1) A is a normal subgroup of E. (2) G = E/A. El- 
ements of G can be viewed as cosets Au{g) in E, where 
u{g) e -E is a representative of g. The choice of rep- 
resentative is arbitrary; that is, we are free to redefine 
u'{g) = a{g)u{g) where a{g) € A. A general element of 
E is of the form au(g), where a G A. This leads us to the 
third and final condition defining an A^-extension: (3) 
uig)a = a<s)u(g)_ 

We have u{gi)u{g2) = ^{gi, g2)ui9i92), where 
'^{91^92) € A is an Ay-factor set again satisfying the 



associativity constraint Eq. (23). At this point it is clear 
that all the structure of factor sets and their equivalence 
classes is identical to the discussion given above, and that 
we can also view these classes as equivalence classes of 
Ay-extensions. 

If G contains no anti-unitary operations, or, more im- 
portantly for the purposes of this paper, if A = Z2, then 
condition (3) above reduces to the statement that A lies 
in the center of E. Such Ay-extensions are called A- 
central extensions. 

For the most part, it is not necessary to use the ter- 
minology of group extensions; we can just as well talk 
about projective representations and factor sets. One 
advantage of the above more abstract discussion is that 
there is no requirement that A be a subgroup of U(l); it 
can be any Abelian group. 

Coming back to projective representations, we see that 
any projective representation belongs to a cohomology 
class. For each class, there are one or more unitarily 
inequivalent irreducible representations. Therefore, clas- 
sifying projective representations by cohomology class is 
coarser than classification by unitary equivalence. 

We now consider a few simple examples to get a feeling 
for the general structure we have been describing. In 
these examples we choose F(l) = 1; this can always be 
done and implies w(l, 1) = a;(g, 1) — a;(l,.g) = 1. When 
doing practical calculations for discrete groups it is often 
convenient to specify the group in terms of generators 
and relations. For instance, G = Z2 is generated by a, 
subject to the relation a? = \. If we consider A = Z2 and 
a projective representation F, the single relation becomes 
[F(a)]~ — a — ±1. Specifying the relation in this way 
defines a factor set uj{1, 1) = Lo{l,a) = Lu{a, 1) = 1, and 
uj{a,a) — a. There are two cohomology classes labeled 
by a, and iJ^(Z2,Z2) = Z2. Each class has two one- 
dimensional irreducible representations: Va=\{a) = ±1, 
and FCT=-i(a) = ±z. 

For any discrete group we can follow this procedure 
of writing down generators and relations. We can write 
the relations so that the right-hand side of each is unity. 
Then, passing to a projective representation F, the right- 
hand side of each relation is replaced by an element of 
A. Let us work out an example to illustrate the proce- 
dure. Suppose G = Z2 X Z2 and A — U(l). We choose 
generators a and h, satisfying the relations 

= 1 (27) 
\? = 1 (28) 
aha-H-^ = 1. (29) 

Passing to a projective representation F, we have 

F(a)' = Oa (30) 
V{hf = at (31) 
F(a)F(6)F(a)-iF(6)-i =aa6, (32) 

where CTq, Oh, Oab e U(l). 

At this point a couple of issues arise. First, the cr's 
are not, in general, in one-to-one correspondence with 
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cohomology classes. That is, there is some redundancy 
that has to be ehminated. Second, some choices of the 
ct's may be inconsistent and not give a legitimate factor 
set. Both these issues arise in this example. We can 
set (Ta , —>■ 1 by redefining the phase of r(a) and T{b). 
Therefore we can write Eq. (32) as 

[r(a)r(&)]2 = aab- (33) 

Multiplying this on the left and right by r(a), we also 
obtain 

[r(6)r(a)]2 = aab- (34) 

Since [r(6)r(a)] = [r(a)r(6)]^^, these two equations are 
only consistent if <Tab = ±1. Therefore we find H^{Z2 x 

Z2,U(1)) =Z2. 

Suppose we consider again G = Z2 x Z2 , but now A — 
Z2. In this case we proceed as above, but cTq, ab, aab G ^2- 
We can no longer eliminate a a and ab, and all choices 
of the ct's are consistent, so we have i?^(Z2 x ^2,^2) — 
Z2 X Z2 X Z2. Notice that the number of classes increased 
upon changing A from U(l) to Z2. Indeed, since every 
Z2-factor set is also a U(l)-factor set, we can group the 
Z2 classes together into U(l) classes: the four Z2 classes 
with aab = 1 belong to the same U(l) class, and similarly 
for the four classes with aab = ^1- 

We can always "coarsen" the Z2 classification in this 
way, grouping Z2 classes together into Ut(1) classes. 
Note the appearance of the T subscript, which is impor- 
tant if G contains anti-unitary operaitons. We denote 
the resulting group of Ut(1) classes by /f^(G,Z2). In 
general, H^{G,Z2) is the subgroup of if^(G,UT(l)) gen- 
erated by all elements of order 2. In the above example 
with G = Z2 X Z2, H'^{G,Z2) = H^{G,\Jt{1)), but this 
is not true in general. We will see that this coarsening 
has an important physical interpretation. 

If G is a continuous group, it is natural that there 
should be some kind of continuity condition on u}{gi, (72 )• 
The naive choice of requiring uj{gi, g2) to be a continuous 
function on the group is not adequate;"*^ for instance, it 
is easily seen that w is discontinuous for the 5 = 1/2 
representation of SO (3). Instead, we believe the correct 
prescription, following Ref. 47, is to require that w ((71,32) 
be a measurable function on G (that is, to classify ex- 
tensions by Borel cohomology). For practical purposes, 
when dealing with continuous groups it is often possible 
to work out the fractionalization classes by simple ele- 
mentary arguments, as illustrated by the discussion of 
Sec. Ill A for G = U(l) and G = S0(3). 

C. General structure, and physical manifestations 
in excited states 

We are now in a position to state the result that frac- 
tionalization classes for each superselection sector are 
given by elements of H'^{G, Z2), where G is the symme- 
try group. We focus on the e-sector only to simplify the 



notation; all statements also hold for m and e sectors. To 
connect with the discussion of the previous two sections, 
we can say that the action of symmetry operations on the 
e-sector states in a region R'^ is given by the projective 
representation F'^, satisfying 

r^(ffi)r^(52) =We(5i,32)F^(5iff2), (35) 

where We ((71,32) G Z2 is a Z2-factor set. 

If a state decomposes into e-sector regions 

Rf, . . . , Wf., then symmetry localization holds, 

Uigm=rig,l)---r%g,km, (36) 

where U{g) is the unitary operator representing g, and 
V^ig, i) is an e-operator on (For a discussion of anti- 
unitary time reversal, see Sec. Ill A.) The notion of e- 
operator was introduced in Sec. II A, and is important 
when the Rf are not connected, which will be the case for 
point group operations as discussed in Sec. Ill D. Physical 
properties are invariant under 

F^(g) ^ A(g)F^(3), (37) 

where X{g) G Z2. This invariance, and the fact that 
'^e(5i,52) S ^2, is a consequence of the fusion rule e x 
e = 1, which also implies k must be even. Due to this 
invariance, the fractionalization class is given by the Z2 
cohomology class of the factor set Wg. 

The fractionalization class is a universal property of a 
Z2 spin liquid phase, so long as symmetry is preserved. 
To sec this, suppose that somehow two e-particles have 
different factor sets uJei and We2, in different classes. 
Then we can fuse them to obtain a 1-particle with factor 
set WeiWe2. But sincc Wei and uie2 are assumed to be in 
different classes, c{uJeiUJe2) 7^ c(l); that is, we have found 
a 1-particle that does not transform in the class of linear 
representations. This is a contradiction, so all e-particles 
must have the same cohomology class. Since cohomology 
classes are discrete, they are then a robust property of a 
phase, so long as symmetry is preserved. 

At this point it is important to ask what type of physi- 
cal information is encoded the fractionalization class, and 
how this information can be extracted. First, mathemat- 
ically, the classification by H^{G, Z2) can be coarsened to 
classification by H^{G, Z2), if we allow for U(l) transfor- 
mations of F"^ (g) [i.e., \{g) E U(l)]. We say that elements 
of H'^{G,Z2) specify the Ut(1) fractionalization class, to 
distinguish it from the Z2 class specified by elements of 
i?^(G, Z2). Physically, U(l) transformations leave mea- 
surable properties invariant in a process during which e- 
particles do not fuse (and are not created in pairs) , except 
for overall phases (and hence eigenvalues) of symmetry 
operators. This can occur, for instance, if we have several 
e-particles that are very far apart and remain far apart 
on some timescale of interest. During such a process the 
number fie of e-particles is a well defined integer con- 
served quantity. The transformation T'^{g) — >■ e^'^r'^{g) 
modifies U{g) — )■ e^'^""'U{g). (The same holds for anti- 
unitary time reversal.) In general, the transformed U{g) 
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is not a symmetry operation, but it is during the process 
of interest (by assumption). Therefore, we can think of 
the Ut(1) fractionahzation class as capturing some prop- 
erties of individual anyons, while the additional informa- 
tion in the TLi class can only obtained when we consider 
fusion of anyons or eigenvalues of symmetry operators. 

It is helpful to illustrate this discussion with an exam- 
ple of G = Z2 X Z2 unitary internal symmetry, which 
also gives some sense of how fractionahzation class in- 
formation may be extracted physically. We take gener- 
ators a and h for Z2 x Z2, satisfying the relations given 
in Eqs. (27)-(29). As discussed in Sec. IIIB, there are 
two Ut(1) classes, depending on whether a and h com- 
mute ((Tab = 1) or anticommute {dah = — !)• Suppose we 
consider energy eigenstates with two localized e-particles 
that are held fixed in space, far enough apart so they do 
not interact with one another. Also suppose that we con- 
sider such states on the sphere, so we do not have to worry 
about the global topological degeneracy. When aab = 1, 
there are four one- dimensional projective irreducible rep- 
resentations. Because the irreducible representations are 
one-dimensional, in the absence of other symmetries, the 
states we consider will be nondegenerate. However, for 
Uab = —1, there is a single two-dimensional irreducible 
representation. This implies that the states we consider 
are fourfold degenerate, because each e-particle has in- 
ternal degrees of freedom described by a two-dimensional 
Hilbert space. 

Using degeneracy of levels works to distinguish the two 
Ut(1) classes, but it does not distinguish the Z2 classes 
within a given Ut(1) class. This makes sense in light 
of the physical interpretation we gave of Ut(1) vs. Z2 
classes; degeneracy of levels has to do with the projective 
irreducible representations associated with individual e- 
particles, but does not involve fusion or eigenvalues of 
symmetry operators. Moreover, given one of the four Z2 
classes within a given Ut(1) class, the other three can be 
realized by making transformations r'^(a) — iT'^{a) and 
y'{h) — )■ ir'^(6); this does not affect the dimensions and 
multiplicities of irreducible representations. 

Extracting the additional Z2 class information is more 
subtle. Suppose we consider the two classes with ab — 
(Jab — 1, and suppose we make the assumption that the 
two e-particles are identical. In the class Ua = 1, in any 
irreducible representation r(a) = ±1, so acting on a state 
with two identical particles we have U{a) — 1. On the 
other hand, in the class with Ua = —1, r(a) — ±i in any 
irreducible representation, so U{a) — —1 on a state of 
two identical e-particles. Subtle information of this kind 
is not captured in the Ut(1) class, as the eigenvalues of 
symmetry operators are involved. Somewhat less obvi- 
ously, fusion is also involved via the implicit assumption 
that the ground state (no e-particles) satisfies U{a) — I; 
in fact, what we are doing is comparing eigenvalues of 
U{a) for the ground state and a state of two identical 
e-particles. This comparison is not well-defined if fusion 
processes are suppressed, because in that case the phase 
of U{a) can be adjusted separately for the two states be- 



ing compared. 



D. Space group symmetry 

Much of the discussion above carries over for general 
space group symmetry, but the notion of symmetry lo- 
calization needs to be modified. This is so because space 
group operations such as reflection and rotation move 
some points by large distances, and can thus move an 
anyon out of the region in which it is localized. We again 
focus on e-particles for concreteness, but all statements 
also hold for m and e particles. 

It is simplest to illustrate the differences from the case 
of translation and internal symmetry by focusing on a 
concrete example. We consider the reflection symme- 
try Px sending x — )■ —x, y y, and satisfying the 
relation = I. As in Sec. Ill A, we consider states 
IV'a) decomposed into two e-sector regions i?^, i — 1,2 
(Fig. 8). Reflection maps these regions to image regions, 
Px '■ Rf ^ Rfi. Symmetry localization can again be ex- 
pressed by writing 

PM = p^a)p^{m<.). (38) 

but now P^{i) is an e-operator on Rf, deflned as the 
union of Rf and i?^,. P^{i) has a single e-string connect- 
ing Rf and Rf, (Fig. 8). Again, the physical interpreta- 
tion is that P^{i) is a "one-particle" symmetry operator. 
This operation can no longer be accomplished entirely 
locally, because the e-particle must be moved from i?f to 
i?^/, hence the presence of the string operator. However, 
once the e-particle is moved, any remaining operations 
can be accomplished locally in Rf and i?^, . 

Just as for the case of translation and internal symme- 
try, we should also consider the effect of any phase factor 
obtained by transforming the string connecting the two 
e-particles. Here, one can follow essentially the same ar- 
gument, illustrated in Fig. 8b, to show that this phase 
factors into a product of phases associated with the in- 
dividual particles. 

Now consider the relation P^ — 1. Arguing as before, 
we have P^il')P^il) = P^{2')Pm = ±h where P|(/) 
is the operator giving the action of Px on the transformed 
e-particle in region R^,. P!j.(i')P^[i) is an e-operator on 
Rt. 

More generally, suppose we consider a group relation 
Si - ■ ■ Sk — 1. For the e-particle in i?f in \ipa), we then 
have Sf ■ ■ ■ = ± 1 , where for simplicity we have sup- 
pressed region labels for the Sf operators. Each of the 
5,f , and thus the product Sf ■ ■ • S^, is an e-operator on 
Rl, defined as the union of Rl, Sk{Ri), Sk-iSk{Ri), and 
so on. Note that we are assuming that the symmetry op- 
erators for one particle commute with those for the other. 
This will be the case if Rl and i?2 do not overlap, which 
we assume. This amounts to assuming that the differ- 
ent particles occupy generic, i.e. non-symmetry related, 
positions. It would be interesting to consider the impli- 
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FIG. 8. (a) Illustration of the initial state \'4>a), and the ac- 
tion of Px on \'4ia)- The vertical dashed line is the reflection 
axis, and the regions Rl, R1, are defined in the main text. 
The e-strings connecting regions for operators Oa and P^{i) 
are shown. The dashed-line O'a string is the image of the Oa 
string under P^. (b) By acting with a closed e-string loop as 
shown, the Oa and Px{i) strings can be eliminated in favor of 
the O'a string. Following the argument depicted in Figs. 7b 
and 7c, the closed loop can be decomposed into smaller re- 
gions with unit eigenvalue acting on the ground state, and 
the effect of transforming the string can be absorbed into the 
definition of Px{i)- In the e-particle case, depending on the 
orientations of the strings in (a), the loop in (b) may need 
to be twisted in some of the regions, to ensure the correct 
orientation of the O'a string. 

cations of relaxing this assumption, but we leave this for 
future work. 

Note that depends on the group relation consid- 
ered. This is unappeaUng, because these regions serve to 
define the e-sector states associated with each particle, in 
which the S'^ symmetry operators act. A solution to this 
is instead to define to be the union of all regions that 
can be obtained as images of Rl under all point group 
operations with some arbitrary fixed center of symmetry. 
We can choose the generators of the symmetry group 



to leave the chosen center of symmetry fixed (or nearly 
fixed). Then, for any relation involving a small number 
of generators, i2f includes i?f as a subset. 

With these modifications, the general structure de- 
scribed in Sec. Ill C continues to hold. In particular, 
fractionalization classes are again given by elements of 
H'^{G, 1^2)1 where G is the full symmetry group includ- 
ing space group operations. 

E. Example: square lattice space group, time 
reversal and spin rotation 

Because it is important for discussing the toric code 
model, and also to make contact with projective symme- 
try group classification, we discuss the example of square 
lattice space group symmetry, combined also with time 
reversal and S0(3) spin rotation. For clarity of notation, 
we focus on the e-sector. 

The symmetry group G is generated by: (1) P^, reflec- 
tion X — )■ —X, y ^ y; (2) P^y, reflection x ^ y; (3) T^, 
translation by one lattice constant along the x-axis; (4) 
time reversal 7"; and (5) Rs{Oh), spin rotation about axis 
n through angle 9. It is convenient to use Ty — P^yT^P^y 
in some of the relations, which is a translation by one lat- 
tice constant along the y-axis. In the e-sector, including 
factors of ±1 to specify a non-trivial factor set, the rela- 
tions are 
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where the cr's are Z2-valued parameters. We also have 

Rl{di^)Q' = a^i?^(M), (40) 

where we can substitute = P-x: Pxy: Px^'T'^- the sec- 
ond set of relations, one might worry that the right-hand 
side can be multiplied by a measurable (but not con- 
tinuous) ±l-valued function /(0), where we must have 
/(O) — 1, since -Rf(O) = 1. However, it can be shown 
that f{9) = 1 for all by assuming i?J(6ln) = e*^-^", and 
solving for 

This is manifestly continuous in 6, and therefore f{9) ~ 
1. 
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The relation Eq. (39k) simply tells us whether e- 
particles carry integer (af^ = 1) or half-odd-integer (cr|j = 
— 1) spin. The other relations are all clearly invariant 
under Z2-valued redefinitions of any of the e-sector gen- 
erators, as in Eq. (37). Moreover, it can be shown that 
all choices of the cr's are consistent. Therefore, we have 
shown that H'^{G,Z2) = ^2^, and there are 2^^ fraction- 
alization classes. If we remove spin rotation symmetry, 
then H^{G, Z2) = Zj", and there are 2^° fractionalization 
classes. 

It is also interesting to work out the Ut(1) fraction- 
alization classes, that is to compute //^(G, Z2). Allow- 
ing U(l) phase redefinitions of the symmetry generators, 
we can choose the phase of P^, P^y and so that 
'^px^'^pxy^^txpx ~^ 1- Upon fixing these parameters, the 
residual phase freedom does not affect any of the other 
relations. The anti-unitary nature of T implies that ad- 
justing the phase of T does not affect any of the rela- 
tions. Finally, crfj is clearly unaffected. Therefore we 
have_i/2(G,Z2) = Z|, or, without spin rotation symme- 
try, H'^{G, Z2) = Zj. The latter result can be extracted 
from Ref. 49, which is a check on the validity of the above 
calculations. 



IV. SYMMETRY CLASSES 

A. General results 

Due to the fusion rule e = e x m, fractionalization 
classes for the three non-trivial anyons cannot be speci- 
fied independently. Instead, knowledge of e and m frac- 
tionalization classes determines the e class. Therefore, 
specifying e and m fractionalization classes specifies a 
symmetry class for a Z2 spin liquid phase. The crucial 
issue, addressed in this section, is to understand how the 
e fractionalization class is determined by the e and m 
classes. 

We now state our results, which we establish in 
Sees. IV B and IV C, where we also provide examples. 
Let the Z2 factor sets associated with the e,m and e 
fractionalization classes be We, i^m and uj^. With only 
translation and internal symmetries, — uif^Um- That 
is. We is given in terms of We and iOm by the i?^(G, Z2) 
group product. In the general case where G includes 
point group operations, then Wf = ajtUJeOJm, where Wt is 
another Z2 factor set depending on the group in a man- 
ner specified below. We refer to the presence of ujt as a 
"twisting" of the 7?^(G,Z2) group product. Physically, 
this twisting is a consequence of the nontrivial braiding 
statistics of e and to, and occurs because products of 
point group operations can braid an e and m bound to- 
gether to form an e, in contrast to translation and internal 
symmetry. 

Whether or not point group symmetry is present, 
the symmetry class can be specified by two elements 
of 77^(G, Z2), one for the e-sector and one for the m- 
sector. Equivalently, we can specify a single element 



of H'^{G,Z2 X Z2). We now discuss the number of dis- 
tinct symmetry classes. The number of fractionalization 
classes is flf — |iJ^(G,Z2)|. Naively wc might say that 
the number of symmetry classes is simply but this 
is not correct, because two classes related by relabeling 
e -f-)- TO are not in fact distinct. This means that ele- 
ments of 77^(G,Z2 X Z2) are not actually in one-to-one 
correspondence with symmetry classes. Taking this into 
account, the number of symmetry classes is 

n,^^{n}-nj) + nf. (42) 

We apply this result to the case of square lattice space 
group, time reversal and spin rotation symmetry, where 
nf = 2", and so flc = 2098176 « 2^1. Or, removing 
spin rotation, Qj = 2^0, so = 524 800 « 2^^. 

B. Translation and internal symmetry 

To relate the e fractionalization class to the e and 
TO classes, we consider e-paricles formed as e-m bound 
states, and compute the fractionalization class of the 
bound state in terms of the classes of its constituents. 
We do this first for the simpler case of translation and 
internal symmetry. 

We consider states jV'a) that can be decomposed into 
two fixed e-sector regions Rj {i = 1,2). Each of these 
regions is further decomposed into an e-sector region (Rf) 
and an m-sector region (i?™)- We assume that 

|Va)=O^OL"lV'o), (43) 

where \ipo) is a ground state. is an e-operator on the 
union of i?f and with an e-string connecting the two 
regions. Similarly, O™ is an w-operator on the union of 
i?™ and i?™. Again, for simplicity but not by necessity, 
we assume \iPq) is a singlet under all symmetry opera- 
tions. 

Now, let Sa (a = 1, 2, 3) be translations or unitary in- 
ternal symmetry operations, satisfying 5*1 5*2 = Sym- 
metry localization for the e-sector regions is expressed by 
writing 

SM=Slil)Sli2M^), (44) 

where S^{i) is supported on The symmetry can be 
further localized to the e and m subregions, that is 

Slit) = Sl{t)ST{t), (45) 

where S^{i) and S^{i) are supported respectively on Rf 
and i?™- The S^{i) are the same operators appearing in 
the localization of Sa if the m-particles are not present, 
and correspondingly for the ^"(i). 
Now, suppose 

St{i:)S-2{t)^oJe{l,2)SI{i) (46) 
^rW^2"W=w™(l,2)^3™«- (47) 
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Then it follows immediately that 

=c.e(l,2K„(l,2)5^(z) =c.,(l,2)5^(z), (48) 
and therefore 



(49) 



the desired result. The same statement holds when we 
consider anti-unitary time reversal symmetry; this is eas- 
ily seen following the discussion of Sec. Ill A. 

We now apply these results to the simple case of trans- 
lation as the lone symmetry. In Sec. Ill A we found 
that there are two fractionalization classes in this case, 
parametrized for the e-scctor by T^T^{T^)-^{T^)-^ = 
'^txtv — ill with corresponding relations for m and e sec- 



^txty<xty At this 



tors. Equation (49) implies crl^ty 

point we might naively conclude there are four symmetry 
classes labeled by ordered pairs (tr^^j^, trj^j^). However, 
the (-1-1,-1) and (—1,-1-1) classes are related by relabel- 
ing e o m, and thus are not actually distinct. Therefore, 
there are three symmetry classes in the case of transla- 
tion symmetry alone. 

We note that there are similarly three symmetry 
classes in the case of SO (3) spin rotation symmetry alone. 
Upon substituting Rl{2TTn) — cr|j for the translation 
symmetry relation, and similarly for m and e sectors, 
the discussion above holds unchanged. 



C. Space group symmetry 

When the symmetry group includes point group oper- 
ations, the result Eq. (49) is modified due to the mutual 
statistics of e and m particles. This leads to a twisting of 
the group product giving the e fractionalization class in 
terms of the e and m classes; in particular, = wtWeOJm, 
where Wt is another Z2 factor set encoding a twisting of 
the i?^(G, Z2) group product. We are interested in de- 
termining the cohomology class of Wt, and thus the e 
fractionalization class. 

We proceed by considering symmetry operations Sa 
(a = 1, . . . , fc), satisfying the group relation Si ■ ■ ■ = 1. 
We suppose that Sf-'-Sl = and S^^ • • • 5^ = cr™, 
and would like to compute S\ - ■ ■ = (j^- We consider 
states as described above in Sec. IV B. Following 

Sec. Ill D we let R\ be the union of images of R\ under 
arbitrary point group operations (with fixed center of 
symmetry), and similarly for i?™ and i?™. We will see 
that CTg = GtOi^Urm whcrc Oi enters as a product over three 
different types of statistical phase factors. Knowledge of 
Ot for enough group relations determines the factor set 



First, we consider a single operation Sa and examine 
the statement of symmetry localization. We again have 



= Sl{\)Sl{2)\i,^) 



(50) 



However, Eq. (45) no longer holds, and is instead modi- 
fied to 



Sl{i) = I¥,(^)5:(^)5^(^) 




(51) 
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FIG. 9. (Color online) Computation of W(i) phase factor for 
P:c reflection symmetry, (a) The initial state \'^a) is depicted 
on the right-hand side, while the vertical dotted line is the 
reflection axis, and the flnal state Px\'>l>a) lies to the left of 
the axis. Regions 7?| and their images under P^, are indicated 
with dotted lines. All particles and strings are shown in black 
in the initial state, and in gray (blue online) in the flnal state, 
e-particles are flUed circles, m-particles are crosses, e-strings 
are solid lines, and m-strings are dashed lines, (b) Depiction 
of the state PJ(l)Pr(l)PJ(2)Pr(2)|V'c.). To compute W(i), 
we simply bring the strings into the flnal state conflguration 
shown in the left-hand side of (a). We flnd W(X) ~ —1, where 
the minus sign arises from crossing the m-string beneath the 
e-string in Rl. We also flnd W{2) — —1, where in this case 
the sign arises from sliding the e-string over the flnal-state m 
particle in Pj'- 



where Wa{i) = ±1 is a statistical phase factor originating 
from anticommutation of e and m strings. The factor 
Wa{i) can be computed, simply by comparing Salipa) 
with S'^(l)5^(l)5'^(2)5^'(2)|V'a). We give an example 
of such a computation in Fig. 9. 
So we have 
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a—1... .,k 



(52) 
(53) 



Here, £{i) = ±1 arises from the fact that sonie of the Sl{i) 
and S^{i) may anti-commute due to crossings of strings. 
In more detail, we observe that S^{i) and Yl^ S^{j) 
define closed loops of e-string and m-string, respectively, 
with pieces of string running among the components of 
Rl and i?™. These two loops are assembled "piece-by- 
piece," reading from right-to-left in the product of sym- 
metry operations in Eq. (52); graphically, strings fur- 
ther to the left in the product can be drawn on top of 
strings further to the right. The factor simply mea- 
sures the Z2 linking number of these two loops; that is 
= (—1)"=, where ric is the number of times m-strings 
need to be crossed below e-strings so that the m-string 
loop lies entirely underneath the e-string loop. 

To compute now we need only act on j-^a) with the 
products of symmetry operations in Eq. (53). We have 



Sl{i)---Sl{i)\i;, 
ST{^)■■■STm, 



{i)<7e 



(54) 
(55) 



where Z'^™'{i) = +1 (—1) if an even (odd) number of 
TO-particles in li^a) are enclosed in the e-string loop de- 
fined by Sf{i) ■ ■ ■ S^{i), with corresponding definition for 
Z"^'^{i), reversing the roles of e and m. Therefore we have 
found 



at = Z-"-ii)Z"'^i)iii)[ Yl Wail) 



(56) 



To illustrate this discussion, we compute ujt for the 
case of square lattice space group symmetry (plus time 
reversal). First we consider the relation = 1. In 
Fig. 10, by considering a convenient state IV'q): we show 
that 



e\2 



(Pi) 



(57) 



that is, CTpj. = 1. In more detail. Fig. 10a illustrates the 
symmetry localization of Px\ipa), showing that W{i) — 
— 1 for i — 1,2. The same is easily seen to be true for 
the symmetry localization of Px{Px\tpa)), so these factors 
cancel in the computation of CTp^.. Figure 10b illustrates 
the closed e and m string loops obtained when comput- 
ing = 1 in terms of the one-particle operators P^ii), 
P^{i'), and so on. These loops do not link, so £{i) = 1. 
Moreover, the e-loops do not enclose any m-particles, 
and vice versa, so Z™'^{i) = Z^"^{i) = 1. Therefore, by 
Eq. (56), CTpj. = 1. It is important to emphasize that we 
are free to choose the e and m strings of the one-particle 
operators P^(i): P^i^), and so on, to run horizontally as 
shown in Fig. 10. The result should not be affected by 
this choice; we have not proved this in general, but have 
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FIG. 10. (Color online) Computation of ap^, for the group 
relation = 1. The graphical notation used here is intro- 
duced in Fig. 9. (a) The reflection axis is the vertical dotted 
line. The state \ipa) is depicted to the right of the axis, with 
e and m particles and (vertical) strings drawn in black. This 
is almost the same state considered in Fig. 9, but with a sim- 
pler arrangement of strings. The state Pxl'i/'a) is shown to the 
left of the axis in gray (blue online). The horizontal strings 
depict the operators P^(i) and P"^{i). From this figure it 
can be seen that W{i) = — 1 in the symmetry localization 
of Px\tpa), and the same is easily seen to hold in the sym- 
metry localization of PxiPxli^a})- (b) Closed loops of string 
obtained from P^{i')P^{i) and Pr(i')PT{i)- Positions of e 
and m particles in |^c«) are shown. Since these loops do not 
link, and, for instance, each e-loop encloses no m particles, 
we have = Z^^'ii) = Z™'=(i) = 1. 



experimented with other conventions and always find u*^ 
to be unaffected. 

Proceeding along the same lines we find at — 1 for all 
other relations, except the relation (PxPxy)^ — 1, which 
we now consider. Noting that P,r/2 = PxPxy is a 7r/2 
rotation, for simplicity we instead consider the equiva- 
lent relation ^ ^- same result with- 
out making this simplification.) Figure 11 illustrates the 
state |'!/'a)j and the e and m strings in P^ /2(*) ^^'^ -^71/2 (*) 
are chosen to run as shown. For subsequent applications 
of P^/2, the corresponding strings are obtained simply 
by rotation of Fig. 11. Again, we emphasize that we are 
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FIG. 11. (Color online) Symmetry localization of n/2- 
rotation R„/2 on a state The center of rotation sym- 

metry is the solid square. The initial state \'4^a) has e and m 
particles arranged along a line extending below the center of 
symmetry, and e-sector regions as shown. The particles 
and strings of the final state are shown in gray (blue online). 
The angled e and m strings are the strings of the i?^/2(*) 
and i?™/2(*) operators. Inspection of this figure shows that 
W{i) = 1 for i = 1,2. 



L"' 




R' 



FIG. 12. Closed e and m string operators obtained upon ex- 
pressing il^/jlV'a) in terms of products of one-particle sym- 
metry operators acting on \ipa}- The locations of e and m 
particles in \ipa} are shown, and the closed strings are labeled 
as shown. 



always free to choose the one-particle symmetry operator 
strings to run in this fashion. We find W{i) = 1 for all 
four rotation operators in the group relation. 

Evaluating i?^/2lV'a), we find the four closed e and m 
string loops shown and labeled in Fig. 12, obtained as 
products of the one-particle symmetry operators. These 
loops do not link, so = 1. We define closed e loops, 



L| = L^Lf^. L| encloses a single e-particle, but this 
does not affect evaluation of L^l^'a)- Therefore, in de- 
termining the Z'^™{i) and Z™'^{i) phase factors arising 
from loops enclosing particles, we can consider the i = I 
e and m particles separately from the i — 2 particles. For 
i = 1, Lf does not enclose any particles, so Z^™(1) — 1. 
On the other hand, L™ encloses the i = 1 e-particle, 
so ^"''=(1) = -1. Similarly, Z"'^(2) = 1, because 
does not enclose the i = 2 e-particle. Finally, we find 
Z^'^(2) — —1. Combining all the statistical phase fac- 
tors, we find o^^^^^y = — 1, that is 



I pep 



xy) 



pxpxy 



^ px'pxy^ pxpxy ' 



(58) 



We have thus found ujt, and shown it is a non-trivial 
factor set. 

It should not be surprising that the ujt twisting appears 
in the i?^/2 ~ ^ group relation, because during the course 
of this relation the constituent e and m particles in each e 
particle are braided around one another. It is interesting, 
and perhaps surprising, that this twisting seems to be 
unavoidable for a discrete rotation symmetry. It could 
thus be valuable to obtain a deeper understanding of the 
LUt twisting. Finally, we note that the lack of twisting in 
the other group relations is to be expected. For these 
relations, given an appropriate choice of \'4'a) there is no 
relative motion of the constituent e and m particles, and 
thus no way for braiding statistics to enter. 



V. GENERAL ABELIAN TOPOLOGICAL 
ORDERS 

Here, we briefly discuss the extension of our symmetry 
classification to general Abelian topological orders, with 
the restriction that the symmetry group G consists only 
of translations and internal symmetries. For simplicity 
we restrict to unitary internal symmetries, and discuss 
inclusion of anti-unitary time reversal at the end of this 
section. 

We find the symmetry classes are labeled by elements 
of H^{G, r2), where r2 is the group of fusion rules. This 
agrees with a result asserted by Kitaev.^^''^^ Not all el- 
ements of H'^{G,T2) describe distinct symmetry classes; 
elements related to others by a relabeling of anyons leav- 
ing the topological structure invariant (e.g., e -k^ ra for 
Z2 topological order) correspond to the same symmetry 
class. 

The fusion group r2 is a finite Abelian group, and is 
thus isomorphic to a product of cyclic groups. Suppose 
there are p cyclic factors and r2 = x x • • • x Zfe^. 
Let ei (i = 1, . . . ,p) be the generators for these factors. 
Physically, the e^ form an elementary set of anyons, from 
which any other type of anyon can be obtained by fusion. 
(In the case of Z2 topological order, k\ — k2 = 2, ei = e, 
and e2 — m.) Therefore, the symmetry class should be 
determined by specifying the fractionalization class for 
each of the e; anyons. 
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Since (cj)'^' = 1, in physical states, anyons of type 
appear in multiples of ki. For translation and internal 
symmetry, the property of symmetry localization is ex- 
pected to hold as above. Therefore, in a straightforward 
extension of the discussion given above for Z2 topologi- 
cal order, the action of symmetry on anyons is a Zfc.- 
central extension of G, and the fractionalization classes 
are given by elements of H^{G,Zk-). Then, because 

H^{G,r2) = H^{G,Zk,) X ••• X H^G,Zk^), (59) 

symmetry classes are labeled by elements of H^{G,r2). 

In the future, it would be interesting to generalize this 
result to the case of full space group symmetry. The 
simplest possibility is that the only modification needed 
for space group symmetry is a twisting of the group 
product giving the fractionalization classes of arbitrary 
anyons in terms of the classes. However, unlike for Z2 
topological order, in general Ci-strings do not commute 
with themselves at crossing points, and this may lead to 
new features in the classification. 

Finally, we discuss inclusion of anti-unitary time re- 
versal. First, if all anyons have ki = 2, then no mod- 
ification of the above discussion is needed to incorpo- 
rate time reversal, because anti-unitary complex conju- 
gation acts trivially on elements of Z2 € U(l). On the 
other hand, if some ki > 2, we might imagine that we 
need to account for non-trivial action of complex con- 
jugation on elements of Zfc. G U(l). This is true, but 
is not sufficient; all such cases are beyond the scope of 
our classification because the assumption that time re- 
versal does not permute different types of anyons is ac- 
tually inconsistent with the topological order. To see 
this, if ki > 2, the e| anyon must have either non-trivial 
self statistics, or non-trivial mutual statistics with some 
other anyon. Otherwise, = 1, a contradiction. Let- 
ting the 9s2 be the self-statistics angle of e|, and 6si the 
same for e^, we have, employing the if-matrix Chern- 
Simons approach, 9s2 — 40si Similarly, letting 0,n2 be 
the mutual statistics angle of ef and some other fixed 
type of anyon, and 9mi the same for and the same 
other fixed anyon, we have 6m2 = '^dmi-^'^ Therefore, 
in order for ef to have non-trivial statistics, we must 
have 9ml O;""; or 9^1 ^ 0, 7r/2, tt, 37r/2. Suppose 
9si 7^ 0, 7r/2, TT, 37r/2. Then acting with time reversal on 
a pair of anyons gives a pair of anyons with new self- 
statistics angle 9'^ — ~9si 7^ 9si mod 27r, which is only 
consistent if time reversal transforms ej into a different 
type of anyon. If instead 9mi 7^ 0, tt, the same argument 
shows that time reversal must transform or the other 
fixed anyon into a different type of anyon. Therefore, 
if time reversal symmetry is present and some ki > 2, 
time reversal must permute the different types of anyons. 
This is an additional motivation to develop a full symme- 
try classification for Abelian topological orders in future 
work, including symmetry classes "beyond fractionaliza- 
tion" where anyons are permuted by symmetry. 



VI. EXPLICIT REALIZATION: TORIC CODE 

In this section we show how these ideas work out ex- 
plicitly in the example of the toric code.^^ Some related 
prior results were obtained in Ref. 21. In particular, we 
will work out fractionalization and symmetry classes for 
the case of square lattice space group symmetry alone. 
It is straightforward to include time reversal in the same 
discussion, but we omit this for brevity. We can under- 
stand the possible fractionalization and symmetry classes 
this case simply by omitting the relations of Sec. HIE 
containing time reversal or spin rotation, and keeping 
the remaining six relations. The flf — 2^ fractionaliza- 
tion classes are elements of 7J^(G', Z2) = Zj, and there 
are Sic — 2080 « 2^^ symmetry classes. Of these, the 
toric code model realizes three, using two tunable pa- 
rameters. 

Guided by the general discussion of Sec. IV, it is also 
possible to explicitly work out the Wt twisting involved 
in relating the e fractionalization class to the e and m 
classes. We have done this, but do not present the results 
here; this essentially amounts to a more cumbersome rep- 
etition of the general constructions of Sec. IV C. 

The toric code model was introduced in Sec. II B. We 
assume that \Ke\ 7^ |^^m|, to avoid extra symmetry 
present in that case. When it is necessary know the size 
of the system, we take an L x L torus with L divisible by 
4. For explicit calculations, we will take the coordinates 
of vertices s of the form r = {x,y) = {m,n) with m,n 
integers, and the coordinates of faces p will take values 
{x,y) — {m + 1/2, n + 1/2). Throughout this section, 
using language from the Z2 gauge theory description of 
Z2 topological order, we shall often refer to e-particles as 
charges and m-particles as fluxes. 

When Ke > 0, the ground state has Ag = 1, and a 
vertex s with = — 1 is an e-particle (a charge exci- 
tation). When Ke < 0, the ground state has As — — 1. 
Viewed in terms of the > ground state, this is a 
background charge of e-particles. This situation often 
arises in theories of Z2 spin liquids and, when it appears 
for a gauge theory Hamiltonian, is referred to as odd Z2 
gauge theory.^'' With Kg <0, the excited e-particles now 
correspond to As = +1 vertices. Identical considerations 
relate m-particles (flux excitations) to the value of Bp. 
The three different symmetry classes that this Hamilto- 
nian accesses are realized for Ke,Km > 0; Ke,Km < 0; 
and signfCe = — signify- The last case corresponds to 
two different choices for the signs of and Km', these 
are not distinct symmetry classes as they are related by 
relabeling e O m. 

A. Wave functions 

It will be helpful to have explicit forms for the wave 
functions. We build the ground state off of a reference 
state that minimizes the flux term of the Hamiltonian. 
Let Se = signKe and Sm = sign Km- For Sm = 1, the 
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FIG. 13. The state |ref(-l)). Dark links carry = -1, 
others have = \. The square at the center of all figures 
identifies the origin of coordinates and of point group opera- 
tions. 



reference state |ref(l)) will have — \ on all bonds, as 
discussed earlier [see Eq. (8)]. For Sm = —1 the reference 
|ref(— 1)) has a^^i = 1 on horizontal links and a^^, = ±1 
on alternating columns of links; that is, ^, = —1 for 
X ~ x' = odd. This puts one link with ^.i = — 1 on 
each plaquette p; see figure 13. The full ground state is 
then 

|V'0(Se, S™)) = n 72 (1 + SeAs)\TefiSm)), (60) 



with 



Asipo = SelV'o), Bplpo = s^lVo) (61) 



for all s,p. We require these states to be invariant under 
space group symmetry. This is manifest for Sm — I since 
the reference state has full symmetry. For Sj„ = —1, we 
can use the identity 



l[{SeA, 

see 



(a) (b) 
FIG. 14. (a) Electric string, (b) Magnetic string. 



(a) 



(b) 



(c) 



FIG. 15. The action of translation Tx by one lattice spacing 
along X on the two-charge state in (a). Figure (b) shows the 
result of moving the quasiparticles with single spin- flips, and 
the shaded area in (c) is swept out by sliding the strings to 
their final position. In (b), the thicker string is doubled; this 
convention will be used in subsequent figures, except where 
noted. 



which define corresponding two-particle states 

|r?, r^) = r (r^, r^)|Vo), K, r^) = /"(r^ , r™)|V'o), 

(64) 



(62) 

depicted in figures 14(a) and 14(b). 



for any set of vertices G- The second product will im- 
plement a space group operation (via spin flips) on the 
reference state for an appropriate choice of G- We require 
G to contain an even number of vertices so that the sign 
factors Se cancel; in particular, this forces us to take L 
divisible by 4, because the appropriate G for the reflec- 
tion Pxy contains {L/2)'^ vertices, which can be chosen 
as the vertices r — (odd, odd). 

With the ground states in hand we can work out the 
excited states. Excitations come in pairs, connected by 
strings: m-strings, connecting two fluxes, consist of spin 
flips cr^; e-strings, connecting charges, consist of a^. We 
choose conventional contours for strings in initial states 
{i.e., those states on which we will act with some sym- 
metry operation). For simplicity, we take each contour 
C(r, r') to consist of two straight segments at most; start- 
ing from the leftmost particle, the contour first goes right, 
then up or down as needed. The initial-state string op- 
erators are 

/^(r?,r^)- n "rr', /'"(rr,r^")= 11 (63) 

C(r,r') C(r,r') 



B. Single- particle symmetry operators 

Now consider how translation acts on a pair of charges, 
say. We have 

r.|r^,r^) = |r5+x,r^+x) 

_ aV2-yi„z z I e e\ 

~ *m "ri,ri+x'^r2,r2+xl'-li^2/ 

= (s5(;<,r,+x) {s'l>Lr.+^) 1^1, r^): (65) 

where we have dropped the superscript e in many places 
to ease the notation. The crucial second equality comes 
from sliding the vertical segment of the string L*^ over the 
possible background flux; see Fig. 15. The result factors 
into single-particle operators, as we argued it should on 
general grounds. Parallel arguments apply for fluxes, and 
for Ty , so we have identified single-particle operators that 
act on single sites r, 

T^ir) = s^<r+y T^{r) = s ^ a^^^^^. (66) 
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FIG. 16. (a) Under Px, charges and fluxes in the • region 
follow contours above the origin, those in the o region go 
below, (b) The conventions for P^y are analogous. 



(a) (b) (c) 

FIG. 17. The action of P^ on a two-e state. 



Here we use the floor function [-J so that all the phases 
are real. As shorthand we write 

T|(r): si, T™(r): , ^^(r) : s^^, T^{r) : , 

(67) 

when we only need the single-particle phases, since the 
necessary factors of and just lie on the contours we 
choose to represent the operators. 

In fact, since we prefer to work with T^, Px, and Pxy as 
generators, we will want to verify that these expressions 
for Ty agree with the relation Ty — P^yT^P^y , which is 
done below. Note that for fluxes, r takes values in the 
dual lattice, i.e. at the centers of faces, and the links in- 
dicated by the subscripts on in these formulae should 
be thought of as links of the dual lattice. Also note that 
the toric code is special in that a quasiparticle is localized 
to a single vertex or plaquette. To connect to the general 
formalism developed earlier, we should write quasipar- 
ticle symmetry operators that act on regions. For the 
purposes of the present discussion, though, we shall just 
take the region R to consist of a single vertex or plaque- 
tte. 

The point group operations Px and Pxy move the quasi- 
particles over greater distances, and so we need to choose 
conventional contours Cp^ (r, PxX) and so on. We depict 
our choices in Fig. 16. The strings always run along the 
boundary of a square centered at the origin; the direc- 
tion is set by the initial position of the particle, which 
separates into two regions, so that the string never cov- 
ers more than 180° of angle. This choice of contour is less 
than obvious for Px] we have chosen it to simplify com- 
putations of the product (PxPxy)'^, which is a product of 
four 7r/2 rotations. 

We define the action of a symmetry on states to trans- 
form the string operators in the natural way, by trans- 
forming the coordinates of the spin operators in the 
string. Note that under point group operations, this will 
not always carry a conventional string to a conventional 
string. 

Now we can work out the phases that accompany our 
choices of strings for point group operations. The guiding 



(a) (b) (c) 

FIG. 18. The action of P^ on a two-m state. 

principle is that a point group operation moves a quasi- 
particle along the boundary of a square centered at the 
origin. The relevant distance is captured by the function 

d(r) =max(|x|,|y|) (68) 

in terms of the coordinates x, y of the particle. The strat- 
egy is just as for translations: (1) Act on a two-particle 
state with the one-particle strings; (2) identify the region 
over which the strings need to slide in order to arrive 
at the transformed state. The calculation can be done 
graphically, although we show a more analytic approach 
as well. The distance d{r) takes half-integer values for 
fluxes; suitable integer-valued functions are 

dx{r) = Ld(r)J -H0o"(r), dxy{v) = [d{r) \ + 9^,y (v) , (69) 

where 0'^{r) [do^{r)] takes the value 1 on the region 
marked o in Fig. 16(a) [16(b)] and otherwise. 

Consider first the action of Px on a state with 
two charges, connected by a conventional string, as in 
Fig. 17(a). Adding the strings of spin flips to move the 
quasiparticles to their reflected positions produces 17(b). 
Then sliding all the strings to their final positions sweeps 
out an area, shown in 17(c). If the area were odd, the 
wave function would pick up a sign. However, the rel- 
evant area will always contain an even number of back- 
ground fluxes, since it will be symmetric about the axis 
X = 0. Therefore, P^ needs no extra phase beyond the 
string of spin flips: 

Pxir) : 1. (70) 

Next consider Px acting on a pair of fluxes, Fig. 18. 
Again the swept-out area has reflection symmetry, but 
may contain an odd number of background charges, lead- 
ing to a phase factor. All background charges contained 
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(a) (b) (c) 

FIG. 19. The action of Pxy on a two-e state. 



(a) 



(b) 



(c) 



FIG. 20. The action of P^y on a two-m state. 



(a) 



(b) 



FIG. 21. (a) The action of P^yT^Pl^y on a single charge e. 
(b) Sliding the string of (a) to obtain a contribution to the 
single-particle phase of Ty . 



(a) 



(b) 



in the swept-out area are doubled (i.e., come in pairs), 
except those lying on the axis of reflection. Therefore, 
one merely needs to count the number of charges along 
the reflection axis inside the swept-out area. This is 
\dx{vi) — dx{v2)\ if both fluxes are in the same half plane 
y ^ 0, and dx{vi) + dx{Y2) if one particle is at y > and 
one is at y < 0, as in 18. The resulting sign therefore 
factors as 



(71) 



The calculation of signs for P^y is analogous. In 
Fig. 19, we see that we need only count the number 
of background fluxes within the swept-out area along 
the line x = y; all the other fluxes are doubled and 
do not contribute a sign. The sign is therefore set by 
d{ri) + d{r2) mod 2, and factors as 

P;^y(r) : sfi--). (72) 
For fluxes, the result is similar, see Fig. 20, and we find 
P™(r) : s^W. (73) 



Recall that we only want to deal with three genera- 
tors, so we demand that Ty = P^yTxP^y on each sector. 
Consider this sequence of operations, for example Fig. 21. 
Because [P^yY = {P^Y = ^^ '^^ can use Pxy rather than 
P~y for the first operation. The phases from the single- 
particle operators are 



PUy + l,x)T'x{v,^)PL{x,v) 



d(y+l,x) X d(x,y) 



„x+l 



\y\ > \x\ 
\y\ < \x\ 



(74) 



FIG. 22. The strings for the group relation P^^. In (a), recall 
that the thicker strings are doubled; that is, they represent 
two strings acting in the same position. 



We can see from Fig. 21(b) that there is an extra factor 
of Sm from sliding the string when |y| > |a;|, so that the 
appropriate phase is 



(75) 



as expected. The marginal case |y| = \x\ gives the same 
result. The same calculation gives the same result on the 
m-sector. 



Symmetry group relations and symmetry 
classes 



1. Direct computations 

Let us work out some of the symmetry relations. Con- 
sider the relation P^y = 1 as it acts on a single flux. 
In most cases this involves putting down a pair of iden- 
tical string operators, which square to 1 trivially, as in 
Fig. 22(a). However, our conventions imply that some- 
times the symmetry operators involved in the relation 
may enclose some background charge, as in 22(b). The 
relation must be constant on a given superselection sec- 
tor, and this example provides a test of this claim. The 
flux in question is located at r = (—2^,2^) in units of 
the lattice spacing. Then we compute 



^xy\^2^ 2' xy\ 



2- 2i^ 



L2ij L25J+1 



(76) 
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FIG. 23. The strings for the group relation (P^Pxy)^- In (a), 
the darker curves are triple strings. In (b), all segments are 
drawn with the same weight. 

where C is the box drawn. This must expression must be 
evaluated on the ground state. Since C contains an odd 
number of background charges it evaluates to s^, so that 
we find 

iP"lf = 1 (77) 

in all cases, as expected. Note that the background 
charge appears twice in this calculation, once explicitly 
and once in the construction of the single particle sign 
factor, and both are important in order to arrive at a 
consistent answer. 

The calculation is essentially identical for both and 
P^y in both the charge and flux sectors. The other point 
group relation, (PxPxy)'^ = 1, which describes fourfold 
rotations, is the most complicated. We have arranged our 
definitions of the quasiparticle symmetry operators so as 
to simplify the computation of this relation. Our defi- 
nitions are such that the relation always gives a square 
contour that encircles the origin once, see Fig. 23 (some 
sections of the contour are traversed three times in gen- 
eral). In the case of a charge, this square encloses an even 
number of fiuxes, while for a flux it encloses an odd num- 
ber of charges. One can work out that the single-particle 
signs always cancel and contribute nothing. Therefore, 
one flnds that 

iP^P^y)'-^ {PTPTyY = S,. (78) 

The other relations involve translations. Consider flrst 
the translation relation T^TyT-^Ty^ = 1. On both 
charges and fluxes, the single-particle signs do not con- 
tribute, 

Tx{x- l,y)Ty{x- l,y- l)r^^(a;,?/- l)Ty^{x,y) : 

sLyJsL^-iJ^Ly-iJsW ^ 1^ (79) 

so the relation just measures the background charge or 
flux enclosed in the elementary loop around which the 
group relation transports the particle in question. That 
is, 

^x^y^x ^y — ^nn ^x ^y ^x ^y — *e- V"^^ 

Next, consider the relation TyPxT~^P~^ = 1 on the 
flux sector, for example. There are cases where the 



FIG. 24. The strings for the group relation 
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corresponding strings do not enclose any vertices, see 
Fig. 24(a), and others in which they do, see Fig. 24(b). 
For the case with no vertices enclosed, the Px segments 
contribute no net sign (they cancel), so one is left with 

the product of the Ty phases, si^^^^~^^ = Se since 
x = n + 1/2 for a flux, for some integer n. In the case 
where vertices are enclosed, the Px contributions do not 
cancel each other, but rather cancel the enclosed flux, so 
that the result is always 

T™P™r™-ip™-i - Se- (81) 

In the case of a charge, the signs from Ty give unity 
because the charges lie at integer positions, and 

T^P^T^-'P^-' = 1. (82) 

Finally, the relation TxPxTxP^^ = 1 is sufficiently simple 
to work through that we do not discuss it here. 
In the end we arrive at the result 
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(83) 

We have thus shown that, depending on Se and s„n the 
toric code model realizes three symmetry classes. (Recall 
again that Sg = — 1, Sm = 1 is related to Se = 1, s„i = — 1 
by relabeling e o to, so these are not distinct classes.) 

One interesting feature of this result is the modifica- 
tion of the rotation relation {PxPxy) in the TO-sector. 
The difference is essentially geometrical; a square with 
vertices (charges) at the boundary contains an even num- 
ber of plaquettes (fiuxes), but a square with faces at the 
boundary contains an odd number of vertices. 

2. Alternate approach to relations in the flux sector 

The computations above are all simpler in the charge 
sector, because we have chosen the point group gener- 
ators to leave a vertex invariant. We could also have 
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centered the point group on a face of the lattice. In- 
deed, there is a simple group automorphism that swaps 
the two; we can replace by ~ T^Px- This obser- 
vation provides a simple way to use the e-sector results 
of the previous section to obtain the relations for the m- 
sector. First, using the known relations for generators 
Px J Pxy 7 Tx I we simply compute the relations for the new 
generators P^, Pxy,Tx. We find 



p2 


— (7px 


p2 

^ X 


^px^txpx 


p2 

^ xy 


^pxy 


p2 

xy 


^pxy 


{PxPxy) 


— ^pxpxy 


(pp Y 

1 x^ xy j 


— ^pxpxy^txty 


rp rji rp—\rTi — 1 

-^x^y-^x ^y 


— O'txty 


rji rp rp — 1 T"" — 1 
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Tx Px 'Px Px 
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Px Px Px Px 


— ^txpx 


-^y^ X-^y ^ X 


^typx 


T P 

^y^ X^y ^ X 


^typx^txty 



(84) 



To use this to find the m-sector fractionalization class, 
we view Px and Pxy as point group operations centered 
on a plaquette rather than a site. Then, Px and Pxy 
are site-centered point group generators. In terms of 
the plaquette-centered generators, calculation of the m- 
sector relations is identical to the calculation of the e- 
sector generators in the previous section, and we find 
O'txty = Se, with all other parameters equal to unity. 
Passing from plaquette to site centered generators using 
the relations above, we recover the results of the previous 
section. 



VII. QUANTUM NUMBERS OF DEGENERATE 
GROUND STATES 

In Z2 topologically ordered phases, it is well known 
that the fourfold-degenerate ground states on a torus 
can have different symmetry quantum numbers. As long 
as the symmetry is preserved, any discrete information 
contained in these quantum numbers is a robust, uni- 
versal property of a Z2 spin liquid phase. Even though 
the ground states belong to the 1-sector, their quantum 
numbers can be partially determined given the symmetry 
class. Here, we do this for the simple case of translation 
symmetry alone. The analysis here can be generalized to 
other symmetry groups; it may be useful to do this in 
future work. 

Consider a finite-size system with periodic boundary 
conditions, and translation symmetry generated by Tx 
and Ty. We do not assume any particular Bravais lat- 
tice, and X and y are just labels for two primitive lattice 
translations. We suppose the system has linear dimen- 
sions {Nx,Ny), where for instance Nx is the number of 
primitive cells in the x-direction. 

The crucial observation is that the generators of the 
loop algebra can be viewed as translation of an anyon 



around a loop. That is, we make the associations 

CI = {T;r- (85) 

Crn = (^T;^)N., (86) 

where fi = x,y. This suggests that, for instance, 

TxClT^' = T,^(T;)^- {T^)-\ (87) 

and so on. 

At this point, we can proceed to consider the three 
distinct symmetry classes. In each case we determine 
the relative crystal momenta of the four ground states 
from the above relations. We have also checked these 
results in the toric code model by direct calculation of 
the ground state quantum numbers. 

In the class a^xty — o'tSty = 1, we have 

T^CIT-' - £1 (88) 
T^£™r-i=/:;", (89) 

independent of {Nx,Ny). Suppose |i/;o) is the ground 
state satisfying >C^|?Ao) = IV'o)- This state must have a 
definite crystal momentum since translations commute 
with A basis for the ground state subspace is 

given by {\^o),C7\^o),C'i'\^o),C7C'^\^o)}, and clearly 
all these states have the same crystal momentum as 
\ipo). So we have determined that in this symmetry 
class, all four ground states have the same crystal mo- 
mentum. While this crystal momentum is not deter- 
mined from the present considerations, the relative crys- 
tal momenta of the ground states are determined (and are 
zero). This symmetry class is realized in the toric code 
for Ke,Km > 0, where it is straightforward to find that 
all four ground states have crystal momentum fc = 0. 

Next, we consider the class a-fxty = 1: ^Ixty = — 1- This 
class is realized in the toric code when Kg < 0, K^^ > 0. 
If {Nx,Ny) = (even, even), then everything proceeds as 
above and the ground states all have the same crystal 
momentum. In the toric code all four ground states have 
crystal momentum fe = 0. 

In the case {Nx,Ny) = (odd, even), we have 

TxC^T-' - (90) 
r,£^T-i - (91) 

while both and Ty commute with £™ and We 
can choose \ipo) as above, but now we see that \ipo) 
and £™|-0o) have the same crystal momentum k, while 
the £™|-i/;o) and £™£JJ'|?/'o) have crystal momentum k + 
(0, 7r). In the toric code we find two ground states with 
crystal momentum zero, and two with (0,7r). The re- 
sults are the same when {Nx,Ny) — (even, odd), except 
of course that the relative crystal momentum between 
pairs of ground states becomes (7r,0). 

For {Nx,Ny) = (odd, odd), Eq. (86) actually implies 

that is, a m-particle translated around the "boundary" 
of the system feels a net 7r-flux inside. This is not a 
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consistent state of affairs on a periodic torus, but it can 
be repaired if a single e-particle is forced into the sys- 
tem. Since this e-particle has no preferred spatial posi- 
tion, we expect the excitation spectrum to be gapless for 
{Nx,Ny) — (odd, odd). This is precisely what occurs in 
the toric code model. 

Finally, we consider the class cr^^ty — ^txty — This 
class is realized in the toric code when K^^K^ < 0. 
To analyze this case, we note that o^^^y = 1. There- 
fore we can simply repeat the analysis above for the 
'^txty = ~'^txty class, but Substituting £^ for In par- 
ticular, we choose jV'o) so that C'^^\ipo) = |V'o)> and |V'o) 
has a definite crystal momentum since £^ commutes with 
translations. Depending on {Nx,Ny), the ground states 
thus have precisely the same relative momenta as in the 
'^txty = ~'^txty class. Morcovcr, we also obtain the same 
results from the toric code model. 

The above analysis points out that ground state quan- 
tum numbers do not completely determine the symmetry 
class, even in the simple case of only transation symme- 
try. A very interesting problem for future work is to 
devise a means to completely determine the symmetry 
class entirely from the ground state wavef unctions. 

VIII. COMPARISON TO PROJECTIVE 
SYMMETRY GROUP CLASSIFICATION 

Here, we compare our symmetry classification with the 
projective symmetry group (PSG) classification of par- 
ton mean- field theories for spin liquids. First, we very 
briefly review PSG classification in the setting where it 
was introduced, namely the fermionic parton approach 
to square lattice S = 1/2 Heisenberg models. (See also 
Ref. 27 for a more extended discussion.) We consider a 
system of 5 = 1/2 spins on the sites of the square lat- 
tice, and assume square lattice space group, time reversal 
and S0(3) spin rotation symmetries. The spin operator 
at site r is written as a bilinear of 5 = 1/2 fermionic 
partons, 

with the local constraint of one fermion per site. Defining 

it is straightforward to show that is invariant under 
local SU(2) gauge transformations 

-> GriJr, (95) 

with Gr e SU(2). 

To proceed, one writes down a mean-field Hamiltonian 
^MFT quadratic in the partons. To describe a spin liquid, 
the mean-field theory should respect the full symmetry 
group. In order to leave "Hmft invariant, symmetries are 



in general accompanied by non-trivial gauge transforma- 
tions. For instance, if S* : r — >■ S'(r) is a space group 
operation, then in general 

S -.^r-^G^^Sir), (96) 

where Gf G SU(2). This is permitted because the phys- 
ical spin operators retain the correct transformation law 
5 : 5'r — > 'S's(r). Such action of symmetry on the partons 
is a projective representation of the symmetry group, and 
this projective representation is referred to as a PSG. 

To fully specify a PSG, it is not enough merely to spec- 
ify the action of the symmetry group on the fermions. 
One must also specify the subgroup of gauge transforma- 
tions leaving Hmft invariant. This subgroup is referred 
to as the invariant gauge group (IGG). We restrict at- 
tention to the case of IGG — Z2, since in this case one 
obtains a Z2 spin liquid upon going beyond mean-field 
theory (see below). Such PSGs are referred to as Z2 
PSGs. The nontrivial IGG transformation is ipr — —ipr, 
which clearly commutes with all symmetry operations. 

PSGs can be classified up to unitary equivalence un- 
der SU(2) gauge transformations, Eq. (95); this provides 
a symmetry classification of mean-field parton Hamilto- 
nians with fixed IGG. This is so because the PSG (and 
IGG) can be determined from any Hmft invariant under 
the symmetry group. Then, keeping the IGG fixed but 
otherwise adding arbitrary symmetry-preserving pertur- 
bations at the mean-field level, the PSG, which is dis- 
crete, remains unchanged. In Ref. 20, Wen found 272 
distinct Z2 PSGs on the square lattice. Actually, there 
are a total of 280 distinct Z2 PSGs.^^ Wen assumed that 
spin rotations are not accompanied by any gauge trans- 
formations. In Ref. 27, together with G. Chen, we showed 
that this assumption can be relaxed, leading to 8 more 
PSGs. 

We note that many of these 280 PSGs do not lead to ef- 
fective low-energy theories for gapped Z2 spin liquids. ^° 
Some PSGs do not actually admit a mean-field Hamil- 
tonian with IGG = Z2; such PSGs are said to be only 
"algebraic PSGs," and not "invariant PSGs."^" For ex- 
ample, we found 8 new algebraic PSGs with Chen in 
Ref. 27, but only 4 of these are invariant PSGs admit- 
ting a mean-field Hamiltonian. Moreover, even among 
invariant PSGs with IGG = Z2, the PSG may require 
the fermions to be gapless at one or more points in the 
Brillouin zone. For instance, this happens for all 4 invari- 
ant PSGs of Ref. 27. In such cases one obtains gapless Z2 
spin liquids, to which our classification does not apply. 

The above discussion does not address the question of 
symmetry classification beyond the mean-field level. In- 
deed, to connect to our classification, we need to go be- 
yond parton mean-field theory. For a mean-field Hamil- 
tonian with IGG = Z2, this can be done by minimally 
coupling the fermions to a dynamical Z2 gauge field (see, 
e.g., Ref. 27 for an example of this procedure). If Hmft 
endows the fermions with a gapped excitation spectrum, 
and if the Z2 gauge field is in its deconfined phase, then 
we obtain a low-energy effective theory for a gapped Z2 
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spin liquid. The mean- field fermions are promoted to e- 
particles, and the m-particles are the gapped fluxes of 
the Z2 gauge field. 

We are now in a position to compare PSG classifica- 
tion with our classification. First, as already discussed in 
Sec. I A, PSG classification is tied to parton formalism, 
and we feel that parton formalism is the wrong language 
with which to classify Z2 spin liquids. Moreover, we know 
of no argument that PSG classification continues to hold 
beyond the mean-field level. Next, it is already clear 
that PSG classification does not provide any information 
about the fractionalization class of the m-sector, and thus 
does not give a complete symmetry classification for that 
reason alone. However, it is interesting to discuss the re- 
lationship between the two classifications for the e-sector. 

From the discussion above, we see that a Z2 PSGs is 
a Z2 central extension of the symmetry group; the group 
A is A = — IGG. The e fractionalization class is sim- 
ply the cohomology class of the PSG. The cohomology 
classification is coarser than PSG classification, because 
two unitarily inequivalent PSGs may belong to the same 
cohomology class. Indeed, we have found instances on 
the square lattice where two unitarily inequivalent PSGs 
belong to the same cohomology class. However, in all 
cases we have found, one of the PSGs in such a pair re- 
quires the fermions to be gapless, and our considerations 
do not apply. While we have not found a case where 
two inequivalent PSGs for gapped Z2 spin liquids are 
equivalent under cohomology classification, we have not 
searched systematically for such examples, and it could 
be interesting to do so. 

The statement that cohomology classification is coarser 
than PSG classification might seem somewhat puzzling, 
since there are 2^^ = 2048 distinct e fractionalization 
classes, while there are only 280 distinct PSGs. This 
occurs because S" = 1/2 fermionic partons are not ca- 
pable of realizing every cohomology class. For instance, 
restricting to point group and time reversal symmetry 
(i.e., symmetries leaving a lattice point fixed), there are 
cohomology classes where the smallest irreducible repre- 
sentation has dimension 4, but the 5 = 1/2 partons only 
provide a two-dimensional on-site Hilbert space. 

In discussing PSG classification, we have focused on 
one particular parton representation of S = 1/2 spin 
models. Other parton representations also exist; for in- 
stance, we could have just as well considered S* = 1/2 
bosonic partons, and discussed the PSG classification 
in that case.^^ For every distinct parton representation, 
the PSG classification needs to be redone. Within the 
framework of PSG classification, it is not obvious how 
to compare PSGs obtained using different parton rep- 
resentations. On the other hand, our classification can 
be applied to effective theories for gapped Z2 spin liq- 
uids obtained from any parton construction. If two such 
effective theories belong to different symmetry classes, 
then they describe different Z2 spin liquid phases. On 
the other hand, if two apparently different such effective 
theories belong to the same symmetry class, they may 



describe the same Z2 spin liquid phase. 



IX. DISCUSSION 

We conclude with a discussion of open issues and fu- 
ture directions. Our approach to symmetry classifica- 
tion can likely be extended to arbitrary Abelian topolog- 
ical orders, including space group symmetry. It would 
also be interesting to consider extension to non-Abelian 
topological order. For two-dimensional Z2 spin liquids, 
the problem of full symmetry classification, where some 
symmetry operations may exchange e and m particles, is 
still open. There is also, of course, the problem of full 
classification of symmetric Z2 spin liquids {i.e., beyond 
symmetry classification). Here, X- matrix Chern-Simons 
approaches may prove useful. ^^''"'^ 

Extending symmetry classification to three dimen- 
sional (d = 3) Z2 spin liquids could be particularly inter- 
esting. First, we point out a connection between our d — 
2 classification and the classification of d = 1 symmetry- 
protected topological (SPT) phases.'^^'^^-ss The d = 1 
SPT phases are essentially classified in terms of projec- 
tive transformations of point objects bound to the ends 
of the d = 1 system {i.e., end states). Quite similarly, our 
d = 2 symmetry classification for Z2 spin liquids is based 
on projective transformations of point objects tied to the 
ends of fluctuating one-dimensional strings, namely the 
anyons. The mathematical consequence of this connec- 
tion is the appearance of the second cohomology group in 
both classifications. Now, in deconfined Z2 gauge theory 
in three dimensions, the topological excitations are point- 
like Z2 electric charges and extended vison loops. A vison 
loop can be viewed as the boundary of a highly fluctu- 
ating two-dimensional membrane (with vanishing surface 
tension). By analogy, we then speculate that there is a 
close connection between the "fractionalization class" of 
a vison loop — assuming it can be defined — and the clas- 
sification of c? = 2 SPT phases. 

Returning to the present classification, we give an ar- 
gument that some of our symmetry classes cannot be 
realized strictly in two dimensions. The argument fol- 
lows Ref. 56, which studied surface theories for three- 
dimensional SPT phases. Suppose G is an internal sym- 
metry, and suppose both the e and m particles have a 
fractionalization class admitting no one-dimensional irre- 
ducible (projective) representations. An example is time 
reversal symmetry with (T"^)^ = (T™)^ = —1. Strictly 
in two dimensions, with only internal symmetries, one 
expects there to be a trivial phase with no topological 
order or spontaneously broken symmetry {e.g., a dimer- 
ized phase). To destroy the topological order, one can 
condense either the e or m particle (e is a fermion and 
thus cannot be condensed), but in this situation this must 
be accompanied by spontaneous breaking of G-symmetry. 
Therefore, there seems to be no way to leave the Z2 spin 
liquid and enter a trivial phase, so we expect this situ- 
ation cannot be realized strictly in two dimensions. On 
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the other hand, Ref. 56 showed that such situations can 
be reaHzed on the surface of a d = 3 SPT phase. This 
discussion estabhshes a connection between our classifi- 
cation and the classification of d = 3 SPT phases, which 
could be interesting to pursue in future work. 

Along similar lines, there is another constraint on sym- 
metry classes among certain (strictly two-dimensional) 
models. Consider a model with translation symmetry 
and S0(3) spin rotation symmetry, with an odd num- 
ber of S* — 1/2 moments per unit cell. For such a 
model a trivial gapped quantum paramagnet with no 
spontaneous symmetry breaking and no topological or- 
der is impossible. This implies, for instance, that both 
e and m particles must have non-trivial fractionalization 
classes. If one of these particles had the trivial fraction- 
alization class, upon condensing it one would obtain a 
trivial quantum paramagnet in contradiction to the the- 
orem of Ref. 57. All this discussion points out that it 
is desirable to obtain a better understanding of which 
symmetry classes can occur in various settings. 

Eventually, we hope our results may lead to the de- 
velopment of tests to distinguish different types of Z2 
spin liquids in numerical studies. There is, of course, the 
connection between symmetry classes and ground state 
quantum numbers discussed in Sec. VII. In cases where 
one has numerical access to the excitation spectrum, one 
could potentially obtain information about the "coars- 
ened" Ut(1) fractionalization classes from multiplicities 
(or, more generally, decomposition into irreducible rep- 
resentations) of nearly degenerate energy levels, as dis- 
cussed in Sec. IIIC for the simple example G = Z2 x Z2. 
More ambitiously, it would be interesting and potentially 
useful to understand how to fully determine the symme- 
try class given only the ground state wavefunction(s). 

While we know of no current candidate materials for a 
gapped Z2 spin liquid, it would nonetheless be interesting 
to devise experimental measurements of symmetry class 
information. Thinking along these lines could lead to 
new experimental tests for fractionalization, which could 
potentially be applicable more broadly, for instance to 
gapless spin liquids. 

Finally, we note that Mesaros and Ran have very re- 
cently proposed a classification of topologically ordered 
phases with on-site symmetry.^^ We also note very recent 
related results of Hung and Wen.^^ It will be interesting 
to understand the relationship between our classification 
and these results. 
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Appendix A: Non-singlet ground states 

When discussing symmetry localization [see, e.g., 
Eq. (9)], we assumed we can find a ground state \ipo) 
that is a singlet under all symmetry operations. Here we 
describe how this assumption can be relaxed; this enables 
us to consider the full four-dimensional ground state sub- 
space on the torus. 

We consider states \ipa) where all excitations, and all 
strings separating anyons, are confined to a large box- 
shaped region R as shown in Fig. 25. We take R to cover 
almost the entire area of the system. For such states, we 
choose the non-contractible loops of C%, £™, and so on, 
to run in the space outside R. Restricting to such states 
allows us to break the Hilbert space into "global sectors" 
associated with the four degenerate ground states in a 
well-defined fashion. Formally, we decompose the Hilbert 
space as a tensor product Hb. ® Hg, where Hu is the 
Hilbert space of excitations contained in R, and Hq is 
the four-dimensional Hilbert space of degenerate ground 
states. If O is supported on R, then in the above Hilbert 
space decomposition we write O — Or^^Ig- On the other 
hand, if O is a loop algebra operator, then O = In 
Og- This means that acting on \ipa) with any operator 
supported on R does not affect the global sector degrees 
of freedom. Conversely, acting with any loop algebra 
operator leaves local properties in R unaffected. 

This discussion motivates a generalized version of the 
symmetry localization assumption. Namely, given a sym- 
metry operation Sa, we assume 

Sa\^a)=Sa{R)S^\^c.)- (Al) 

Here, Sa{R) is a unitary operator supported on R. The 
operator is a unitary linear combination of prod- 
ucts of loop algebra generators. Any unitary transfor- 
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mation on the four-dimensional ground state subspace 
can be written as such a hnear combination, so can 
be thought of as a general unitary transformation among 
the global sectors. Clearly we have [Sa{R), S^] — 0. 



At this point we apply symmetry localization as dis- 
cussed previously to the operator Sa{R)- To illustrate 
this with a concrete example, suppose that two localized, 
isolated e-particles are contained within R, in regions Rf 



(i = 1,2). Then we write 

5a|Va) = Sa{R)S^\ij^ = 5a (l)5a (2) | V^„) . (A2) 

Considering a group relation 5i • • • 5*;; = 1, we note that 
we must have 



(A3) 



because Sa — on the ground state subspace, and the 
ground states of course do not transform projectively. 
Therefore the operators drop out in the group rela- 
tions, and discussion of fractionalization classes proceeds 
exactly as in Sec. III. 
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